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ABSTRACT 



We present the second catalog of high-energy 7-ray sources detected by 
the Large Area Telescope (LAT), the primary science instrument on the Fermi 
Gamma-ray Space Telescope (Fermi), derived from data taken during the first 
24 months of the science phase of the mission, which began on 2008 August 4. 
Source detection is based on the average flux over the 24-month period. The 
Second Fermi-LAT catalog (2FGL) includes source location regions, defined in 
terms of elliptical fits to the 95% confidence regions and spectral fits in terms of 
power-law, exponentially cutoff power-law, or log-normal forms. Also included 
are flux measurements in 5 energy bands and light curves on monthly intervals 
for each source. Twelve sources in the catalog are modeled as spatially extended. 
We provide a detailed comparison of the results from this catalog with those from 
the first Fermi-LAT catalog (1FGL). Although the diffuse Galactic and isotropic 
models used in the 2FGL analysis are improved compared to the 1FGL catalog, 
we attach caution flags to 162 of the sources to indicate possible confusion with 
residual imperfections in the diffuse model. The 2FGL catalog contains 1873 
sources detected and characterized in the 100 MeV to 100 GeV range of which 
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we consider 127 as being firmly identified and 1171 as being reliably associated 
with counterparts of known or likely 7-ray-producing source classes. 

Subject headings: catalogs gamma rays: general; PACS: 95.85.Pw, 98.70.Rz 



1. Introduction 

This paper presents a catalog of high-energy 7-ray sources detected in the first two years 
of the Fermi Gamma-ray Space Telescope mis sion by the Large A rea Telescope (LAT). It is 



the succ essor to the LAT B right Source List (lAbdo et al.ll2009dl ) and the First Fermi LAT 



(1FGL, lAbdo et al.ll2010gl ) catalog, which were based on 3 months and 11 months of flight 
data, respectively. The new catalog represents the deepest-ever catalog in the 100 MeV - 
100 GeV energy range and includes a number of analysis refinements. 

Some important improvements compared to the 1FGL catalog are: 

1. The 2FGL catalog is based on data from 24 months of observations. 

2. The data and Instrument Response Functions (IRFs) use the newer Pass 7 event se- 
lections, rather than the Pass 6 event selections used previously. 

3. This catalog employs a new, higher-resolution model of the diffuse Galactic and isotropic 
emissions. 

4. Spatially extended sources and sources with spectra other than power laws are incor- 
porated into the analysis. 

5. The source association process has been refined and expanded. 

Owing to the nearly continuous all-sky survey observing mode and large field of view 
of the LAT, the catalog covers the entire sky with little observational bias. The sensitivity 
is not uniform, due to the large range of brightness of the foreground diffuse Galactic 7-ray 
emission. In addition, because the point-spread function (PSF) and effective area of the LAT 
depend on energy, the sensitivity limit depends markedly on the intrinsic source spectrum. 

As has been established with the 1FGL catalog, a number of source populations are 
known to be present in the data. For individual sources, associations with objects in other 
astronomical catalogs are evaluated quantitatively. 
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In Section 2 we describe the LAT and the models for the diffuse backgrounds, celestial 
and instrumental. Section 3 describes how the catalog is constructed, with emphasis on what 
has changed since the analysis for the 1FGL catalog. The 2FGL catalog itself is presented 
in Section 4, along with a comparison to the 1FGL catalog. We discuss associations and 
identifications in Section 5. After the conclusions in Section 6 we provide appendices with 
technical details of the analysis and of the format of the electronic version of the 2FGL 
catalog. 

2. Instrument & Background 
2.1. The data 



The LAT is a 7-ray detector designed to distinguish 7-rays in the energy range 20 MeV 
to more than 300 GeV from the intense background of energetic charged particles found 
in the 565 km altitude orbit of the Fermi satellite. For each 7-ray, the LAT measures its 
arrival time, direction, and energy. The effective collecting area is ~6500 cm 2 at 1 GeV 
(for the Pass 7 event selection used here; see below), the field of view is quite large (>2 sr), 
and the observing efficiency is very high, limited primarily by interruptions of data taking 
during passage of Fermi through the South Atlantic Anomaly (~13%) and trigger dead 
time fraction (~9%). The per-photon angular resolution is strongly dependent on energy; 
the 68% containment radius is about 0?8 at 1 GeV (averaged over the acceptance of the 
LAT) and varies with energy approximately as E~ - 8 , asymptoting at ~0?2 at high energies. 
The tracking section of the LAT has 36 layers of silicon strip detectors to record the tracks 
of charged particles, interleaved with 16 layers of tungsten foil (12 thin layers, 0.03 radiation 
length, at the top or Front of the instrument, followed by 4 thick layers, 0.18 radiation length, 
in the Back section) to promote 7-ray pair conversion. Beneath the tracker is a calorimeter 
comprised of an 8-layer array of Csl crystals (1.08 radiation length per layer) to determine 
the 7-ray energy. The tracker is surrounded by segmented charged-particle anticoincidence 
detectors (plastic scintillators with photomultiplier tubes) to reject cosmic-ray background 



Atwood et a 
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J2012b) 



(200 9|) and the in -flight calibration of the LAT is described in I Abdo et al. 



The data analyzed here for the 2FGL catalog were taken during the period 2008 August 
4 (15:43 UTC) - 2010 August 1 (01:17 UTC). During most of this time Fermi was operated 
in sky-scanning survey mode (viewing direction rocking north and south of the zenith on 
alternate orbits). Time intervals flagged as 'bad' (a very small fraction) were excluded. Fur- 
thermore, a few minutes were excised around four bright GRBs (GRB 080916C: 243216749- 
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243217979, GRB 090510: 263607771-263625987, GRB 090902B: 273582299-273586600, GRB 090926A: 
275631598-275632048 in order to avoid having these bright transients distort the analysis 
of the more persistent catalog sources near these directional). We are preparing a separate 
catalog of LAT GRBs. 

Previous analysis of the Fermi LAT data relied on criteria for selecting probable 7-ray 
events from all the instrument triggers as determined before launch or modified versions of 
these selections (called Pass 6_V3 Diffus^§). Experience with the data allowed us to develop 
an improved event selection process with lower instrumental background at energies above 10 
GeV and higher effective area at energies below 200 MeV. These Pass 7_V6__ (P7_V6) Source 
class event selections are accom panied by a corresponding revised set of Instrument Response 



Functions (lAbdo et al.ll2012bl ). including an energy-dependent PSF calibrated using known 
celestial point sources. The model for the diffuse gamma-ray background was fit using P7_V6 
Clean event selections and IRFs (see § 12 .2 j) . The Clean event selection has lower residual 
background intensity than P7_V6 Source at the cost of decreased effective area, a tradeoff 
that is worthwhile for studies of diffuse 7-ray emission. The IRFs tabulate the effective area, 
PSF, and energy dispersions as functions of energy and inclination angle with respect to the 
LAT z-axis. The IRFs are also tabulated as a function of the location of the 7-ray conversion 
in the LAT; Front conversions occur in the top 12 tracking layers. The tungsten foils are 
thinnest in this region and the PSF is narrower than for the Back section, which has 4 layers 
of relatively thick conversion foils. The 2FGL catalog is therefore derived from a new data 
set rather than simply an extension of the 1FGL data set. 

During the 1FGL time interval (up to 2009 July 4) the standard rocking angle for 
survey-mode observations was 35°. During much of 2009 July and August it was set to 39°. 
Then on 2009 September 2 the standard rocking angle was increased to 50° in order to lower 
the temperature of the spacecraft batteries and thus extend their lifetime. Time intervals 
during which the rocking angle of the LAT was greater than 52° were excluded. The more- 
conservative 1FGL limit of 43° had to be raised to accommodate the larger standard rocking 
angle. 

For the 2FGL analysis we apply a more conservative cut on the zenith angles of the 7- 
rays, 100° instead of the 105° used for the 1FGL catalog. This compensates for the increased 
contamination from atmospheric 7-rays from the earth's limb due to the larger rocking angle. 
Another motivation for the tighter cut is that the new Pass 7 event selections used for the 



1 These are Mission Elapsed Times, denned as seconds since 00:00:00 UTC on 2001 January 1. 
2 http://www. slac.stanford.edu/exp/glast/groups/canda/archive/pass6v3/lat_Performance. htm 
3 http: / /www. slac.stanford.edu/exp / glast / groups / canda/archive /pass7v6 /lat_Performance.htm 
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2FGL analysis have much greater effective area at low energies than those used for the 
1FGL analysis. Because the point-spread function broadens with decreasing energy, a more 
conservative limit on zenith angle is warranted in any case. 

The energy flux map of Figure [1] summarizes the data set used for this analysis. The 
corresponding exposure is relatively uniform, owing to the large field-of-view and the rocking- 
scanning pattern of the sky survey. With the new rocking angle set to 50° the exposure is 
minimum at the celestial equator, maximum at the North celestial pole and the contrast 
(maximum t o minimum exposu re ratio) is 1.75 (Fig. [5]). The exposure with rocking angle 



35° (Fig. 2 of lAbdo et al.ll2009dl ) was least at the South celestial pole, with a contrast of 1.33. 
The North/South asymmetry is due to loss of exposure during passages of Fermi through 
the South Atlantic Anomaly. Figure [3] shows that the original rocking scheme resulted in 
a very uniform exposure over the sky. The new rocking scheme is less uniform, although 
it still covers the entire sky to an adequate depth. The exposure map for 2FGL is about 
halfway between the 35° and 50° maps. It peaks toward the North celestial pole and is rather 
uniform over the South celestial hemisphere, with a contrast of 1.37. Note that the average 
etendue of the telescope is only slightly reduced, from 1.51 m 2 sr (at 1 GeV) in the first 11 
months to 1.43 m 2 sr over the last 11 months. The reduction is due to the part of the field 
of view rejected by the newer zenith angle selection. 



2.2. Model for the Diffuse Gamma-Ray Background 

The 7-ray emission produced by the Galaxy originating from the interaction of cosmic- 
ray electrons and protons with interstellar nucleons and photons is modeled with the same 
method as for the 1FGL catalog. We fit a linear combination of gas column densities, an 
Inverse Compton (IC) intensity map, and isotropic intensity to the LAT data using the 
P7_V6 Clean data set. To account for the non-uniform cosmic-ray flux in the Galaxy, the 
gas column densities are distributed within galactocentric annuli. More details on the various 
radio and infrared surveys used to generate the maps for the different annuli are given at the 
Web site of the Fermi Science Support Centeifl Inverse Compton 7-rays from cosmic-ray 
electrons interactin g on optical, infrared and CMB photons are modeled with GALPROP 



( jStrong et al.l 120071 ) . In each energy band, the gas emissivities and IC normalization were 



left free to vary. 

For this study we have improved the modeling of the diffuse emission in several ways. 
With more than twice the 7-ray statistics we were better able to discriminate between the 



4 http: / /fermi. gsfc.nasa.gov/ssc / data/access /lat /BackgroundModels.html 
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Fig. 1. — Sky map of the energy flux derived from the LAT data for the time range analyzed 
in this paper, Aitoff projection in Galactic coordinates. The image shows 7-ray energy flux 
for energies between 100 MeV and 10 GeV, in units of 10 -7 erg cm -2 s _1 sr^ 1 . 




0.42 0.84 1.3 1.7 2.1 2.5 2.9 3.4 3.8 4.2 



Fig. 2. — Exposure of the LAT for the period September 2009 to July 2010 when the rocking 
angle was 50°, Aitoff projection in Galactic coordinates. The units are equivalent on-axis 
exposure at 1 GeV in Ms. 
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Fig. 3. — Distribution of the equivalent on-axis exposure of the LAT at 1 GeV. The curves 
show the area of the sky exposed at that depth. The dashed curve is for the first 11 months 
(1FGL: 2008 August to 2009 June) when the rocking angle was 35°, and the full curve is for 
the period when the rocking angle was 50° (2009 September to 2010 July, also 11 months). 
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template maps described above and we were also able to increase the number of energy 
bins from 10 to 14, spanning 63 MeV to 40 GeV. Below 63 MeV, the combined effects of a 
low effective area and increased earth limb contamination owing to the increased breadth of 
the PSF prevent study of the diffuse emission. Above 40 GeV the statistics are too low to 
discriminate between the large number of templates that comprise the model. The quality 
of the determination of the linear coefficients (interpreted as the 7-ray emissivities for the 
gas) was also improved at high energies by using the P7_V6 Clean data set, which has lower 
residual charged particle backgrounds at high energies than P6_V3 Diffuse. For energies 
below 63 MeV or above 40 GeV the diffuse emission model was derived by extrapolating the 
measured emissivities according to a fit of the emissivities in terms of bremsstrahlung and 
pion decay components. 

The spatial resolution of the model was improv ed from 0?5 to 0?125, which is the sam- 



pling of most of the CO survey (IDame et al.l l200ll ). The higher resolution in the fitting 



procedure helps discriminate H2, H I, dark gas, and smoother distributions like inverse 
Compton. For the actual fitting, for computational considerations we sampled the maps 
with 0?25 resolution to derive the emissivities and used the full resolution to reconstruct the 
model from the deduced emissivities. The final resolution of the model is then 0?125. Given 
sufficient statistics this is crucial to discriminate point-like sources and molecular clouds at 
the PSF scale. 

This procedure revealed regions with photon excesses not correlated with gas or tem- 
plates defined by observations at other wavelengths. We found what appear to be two distinct 
origins for the excesses, depending on energy. For both cases we introduce ad hoc 'patches' 
in the diffuse emission model to account for their contributions. The patches are regions of 
spatially uniform intensity whose shapes reproduce the shape of the excesses. The inten- 

Table 1. Additional Components in the Diffuse Emission Model 



Designation 


Center 


Approx. dim. 




Fraction of total 


Fraction of intensity 




(l,b) 


(Z x b) 




intensity 


within patch 


First quadrant and inner 


25°, 0° 


40° x 30° 


1.9% 


1.0% 


13.4% 


Fourth quadrant 


-35, 9 


40 x 30 


2.4 


0.3 


3.8 


Lobe North 


0, 25 


50 x 40 


3.9 


0.1 


6.9 


Lobe South 


0,-30 


50 x 40 


3.7 


0.4 


14.1 



Note. — Description of the additional components added in the Galactic diffuse model. The centers and 
extents are in Galactic coordinates. The extents are approximate because the shapes are irregular, f! is the 
solid angle. To evaluate the fractional intensities, we integrated the intensity above 130 MeV for the First 
and Fourth quadrant patches and above 1.6 GeV for the Lobes patches. 
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sity of the emission associated with each patch is fitted for each energy band together with 
the other templates. The shapes of the patches were chosen to approximately encompass 
regions with an excess of photons of the order of 20% compared to the model outside the 
Galactic plane. Two of the regions have a hard spectrum and are lobe-shaped north and 
south from the Galactic center. This emission was also observed and studied in detail by 



Su et all ( 120101 ). Table [T] summarizes the patches and their contributions to the model. Im- 
ages showing the locations and extents of the patches are available from the Fermi Science 
Support Center at the URL cited above. These regions do not correspond to fluctuations 
in the diffuse emission model. We do not see large regions where the model exceeds the 
observed intensity, and we did not need to use 'negative' patches. Four main regions were 
identified in the first and fourth quadrants, and north and south of the Galactic center. We 
added an extra inner patch to the first quadrant region where the intensity was greater than 
in the rest of the patch. The spectra of the patches were determined in the same way as for 
the other templates by extracting their intensities from fits in each energy bin. 

At lower energies, below a few GeV, an excess of photons seems to be associated with 
the giant radio loop Loop I. The North Polar Spur is clearly visible i n the LAT data and can 



be roughly modeled with the 408 MHz radio map of lHaslam et all (119811 ) as well as a large 
rounded shape filling the Loop. At low energies distinguishing between 7-rays originating 
from Loop I and from larger distances is very difficult near the Galactic plane. It is possible 
that the scaling of the model map for the Galactic inverse Compton emission as well as the 
fitted emissivities of inner Galaxy gas rings are artificially increased in the fitting procedure 
to account for 7-rays produced locally. While keeping the overall residual fairly flat, this 
may bias the diffuse emission spectrum and derived spectra and significances of faint sources 
in a large region of about 100° wide in longitude and 30° in latitude centered in the Galactic 
center. Independent of this effect, other regions are probably inadequately modeled, for 
example the Cygnus region, the Carina tangent, and the Orion molecular cloud; see §3.9. 

The spatial grid of the model now has a bin centered at latitude zero. Previously the 
Galactic ridge was split between two bins with the consequence of flattening the modeled 
ridge and possibly inducing the detection of spurious sources in the Galactic ridge. 

We also created a template for the emission from the earth limb that is not completely 
removed from the P7_V6 Source and Clean data sets at energies below 200 MeV. These are 
7-rays that are in the broad tails of the PSF and so pass the selection cut on zenith angle (see 
§ 12. II) . For the template we used the residuals in the 50-68 MeV energy range and assumed 
that the spatial shape is independent of energy. The very soft spectrum was derived by 
adding this template to the model. The template is specific to the data set analyzed here 
because the residual earth limb emission depends on the orientation of the LAT. 
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The isotropic component was derived for the P7_V6 Source data set by fitting the data 
for the whole sky using the Galactic diffuse emission modeled as above. By construction the 
isotropic component includes the contribution of residual (misclassified) cosmic rays for the 
P7_V6 Source event analysis class. Treating the residual charged particles as effectively an 
isotropic component of the 7-ray sky brightness rests on the assumption that the acceptance 
for residual cosmic rays behaves similarly as for 7-rays; in particular we assume that the 
relative contributions of the Front and Back events to the isotropic intensity are according to 
their relative effective areas. This approximation is necessary in the gtlike analysis described 
in § 13.21 The actual residual background rates for Front and Back events do not in fact 
scale precisely with the (7-ray) effective areas, with the most notable difference being in the 
low energy range <400 MeV for which the background 'leakage' in the Back section of the 
tracker is appreciably greater than for the Front section. This has the effect of decreasing the 
flux measurements at low energies (below ~200 MeV) and hardening the spectra, with the 
greatest effects for low-significance, soft sources. On average the spectral indices for power- 
law spectral fits are hardened by less than half of the typical uncertainty in the measured 
spectral index. 

The models for the Galactic diffuse emission and the isotropic background spectrum, 
along with more detailed descriptions of their derivation, are available from the Fermi 
Science Support Center. 

3. Construction of the Catalog 

The procedure used to construct the 2FGL catalog has a number of improvements 
relative to what was done for the 1FGL catalog. In this section we review the procedure, 
with an emphasis on what is being done differently. 

As for the 1FGL catalog, the basic analysis steps are source detection, localization 
(position refinement), and significance estimation. Once the final source list was determined, 
by applying a significance threshold, we evaluated the flux in 5 bands and the flux history 
(light curve of the integrated flux) for each source. 

Also as for the 1FGL analysis, the source detection step was applied only to the data 
from the full 24- month time interval of the data set. We did not search for transient sources 
that may have been bright for only a small fraction of the 2-year interval. See § 15.2. Ill for a 
discussion of transient LAT sources reported in Astronomer's Telegrams. Analysis of 2FGL 
catalog source variability is found in § 13.61 

The 2FGL catalog is primarily a catalog of point (spatially unresolved) sources detected 
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by the LAT in the 24-month interval. As discussed below, the analysis and catalog also 
include a number of LAT sources that are known to be spatially extended. These sources 
are defined specially in the analysis (see § 13 .4j) but are considered members of the 2FGL 
catalog. 



3.1. Detection and Localization 

Detection of point sources involves iterating through three steps: (1) identification of 
potential point sources, denoted as 'seeds', that have not already been selected in a previ- 
ous iteration; (2) a full all-sky optimization of a model of the 7-ray sky (diffuse emission 
plus sources) including the new seeds to refine their estimated positions and evaluate their 
significances; (3) creation of a 'residual Test Statistic (TS) map'. The TS is evaluated as 
TS = 2 (log £ (source) — log £(nosource)), where £ represents the likelihood of the data given 
the model with or without a source present at a given position on the sky. In each case the 
likelihood i s assumed to have b een maximized with respect to the adjustable parameters of 



the model flMattox et al.Ml996h 



We performed this analysis using the pointlike analysis system, for which the data are 
partitioned by whether the conversion occurred in the Front or Back sections of the tracker 
and binned in energy with four bins per decade from 100 MeV to 316 GeV. For each such 



parti tion, or band, the 7-rays are partitioned according to their HEALPix (IGorski et al. 



20051 ) indices, with the nside parameter chosen such that the angular size of the partition is 
small compared with the PSF for that energy and conversion position. Detailed simulations, 
analytic studies, and adjustments of the bin size have shown that this does not lose precision 
compared with a fully unbinned procedure. 

We discuss each step of the iteration in turn. 



3.1.1. Determination of seeds 
We started with an initial model comprised of the 1FGL catalog of sourc es to which we 



Damiani et al. 


1997; 


(Campana et al 


2008) 



19971; ICiprini et al.l 120071 ) , and a minimal spanning tree-based algorithm 



2008) as described in 1FGL. For the 2FGL catalog analysis, we also included 



in the model 12 spatially-extended sources that have been detected by the LAT; see § 13.41 In 
subsequent iterations, seeds may be added by examination of the residual TS map, described 
below. Since source detection is an integral part of the iteration procedure, the efficiency of 
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the initial seed-determination procedures is not critical. 



3.1.2. All-sky optimization 

We define 1728 circular regions centered on points defined by a HEALPix tessellation 
with nside = 12. All 7-ray data within a 5° radius of each of the points are fit to a model 
including the diffuse components described in § !2.2l and all seeds within a radius of 10°. Each 
region was optimized independently. The parameters included the normalization of each dif- 
fuse component and the spectral parameters of the point sources lying inside the boundaries 
of the HEALPix pixel that defined the region. Since neighboring regions are coupled, shar- 
ing data and sources, we repeated this step until the likelihoods were jointly optimized. For 
some regions along the Galactic plane, convergence required up to 10 iterations. 



For point sources identified as pulsars by LAT phase analysis ( lAbdo et al.ll2010tU2012d ). 
the spectra were fit to a power law with an exponential cutoff; others were fitted to either 
a simple power law, or log- normal (also called log-parabola); the latter was used if it sub- 
stantially improved the overall likelihood. These functio ns are described in § 13.31 Each seed 



was characterized by two versions of the likelihood TS (IMattox et al.lll996f ): one measuring 
the spectral-shape independent measure from independent fits of the fluxes in each energy 
band (TSband), and another which is the result of a fit to the spectral model, {TS mo d e i)- The 
former always will be larger than the latter: the difference is used to decide to switch from 
a power law to a log parabola spectral shape. Seeds with TSband < 10 are eliminated from 
further analysis. The rest are retained in the model for the pointlike optimization. After 
the optimization was complete, those with TS mo d e i > 10 were passed on to the gtlike step 
described below, with the pointlike fit as a starting point. 



3.1.3. Residual TS map 

After the analysis in the previous step converged, we performed a special analysis of 
the full sky to search for missing point sources. A HEALPix tessellation with nside = 512 
is used to define 3.1M points on a 0?1 grid. For each point, we added a new point source 
with a power law spectrum and fixed spectral index 2.0, to the model, and the likelihood 
was maximized as a function only of its flux. The resulting array of values of TS is plotted 
as a sky map. 

Clusters were defined by proximity: a cluster is the set of all pixels that occupy adjacent 
positions. The analysis generated a list of all clusters of such pixels with TS > 10 on the 
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map, used as seeds to be added for the next iteration of the all-sky analysis. We estimated 
the position of a presumed source from the centroid of the pixels, weighted by TS; this 
position was refined later if the seed survived the full analysis. Adding seeds from the map 
was done automatically for Galactic latitudes above 5°; along the Galactic plane the data 
are not always well represented by either point sources or the model for diffuse Galactic 
emission, and we introduced new point sources only if they appeared to be well isolated 
under visual inspection. Figure H] shows the final such map for a region along the Galactic 
ridge. 

In total, 3499 candidate sources were passed to the significance and thresholding step 
of the analysis. 

3.1.4- Localization 

The processing that created the residual TS map used for source detection also per- 
formed local optimizations of the likelihood with respect to the position of each point source, 
using the spectral-shape independent definition of the likelihood, TSband, described above, 
with the rest of the model fixed. The positional uncertainty for each source was estimated by 
examining the shape of the log likelihood function, fitting the distribution to the expected 
quadratic form in the angular deviations from the best fit position. A measure of the quality 
of this fit is the mean square deviation of the log likelihood with respect to the fit on a 
circle of radius corresponding to two standard deviations. For the catalog we tabulated the 
elliptical parameters including the fit position and the fit quality. As in the case of the 1FGL 
catalog, we made two empirical corrections based on comparison with the known locations 
of high- confidence associated sources: multiplied by a 1.1 scale factor, and added 0?005 in 
quadrature to the 95% ellipse axes. This latter is comparable to the spacecraft alignment 
precision requirement of 10". 

We searched for systematic biases in source positions, using comparisons with counter- 
part positions (§ E]). Two cases were considered: (1) sources near the Galactic plane, the 
positions of which might have been suspected to be biased by the strong gradient of the 
intensity of the Galactic diffuse emission, and (2) weak sources near much stronger ones. We 
did not find significant biases in either case. In addition, in Appendix [A] we show that the 
sizes of the localization regions for weak sources are consistent with expectations, as is the 
weak dependence on the source spectrum. 
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3.2. Significance and Thresholding 

To evaluate the fluxes and spectral parameters, and also significances, for the catalog 
we use the standard LAT analysis tool gtlike and associated LAT Science Tools^ (version 
v9r23p0). The localization procedure (§ 13. 1 .4j) provides spectra and significances as well, 
but we do not have as much experience with it so we prefer relying on the standard tools 
whenever possible. This stage of the analysis is similar in principle to what was done for 



the 1FGL catalog ( lAbdo et al.l l2010gT ) . It splits the sky into Regions of Interest (Rol) in 
order to make the log£ (where £ is the likelihood function) maximization tractable, varying 
typically half a dozen sources near the center of the Rol at the same time. (There were 933 
Rols for 2FGL.) This requires an iterative scheme in order to inject the spectra of all sources 
in the outer parts of the Rol. It uses the same energy range (100 MeV to 100 GeV) and 
adjusts the source spectra with positions fixed to the result of § 13.1.41 The same parameters 
are used to refit the diffuse emission model (described in § 12. 2j) to each Rol: normalization 
and small corrective slope of the Galactic component and normalization of the isotropic 
component. We define the Test Statistic TS = 2Alog£ for quantifying how significantly 
a source emerges from the background. The iteration scheme was also identical, as well 
as the threshold at TS > 25 applied to all sources, corresponding to a significance of just 
over 4 a eval uated from the y 2 d istribution with 4 degrees of freedom (position and spectral 



parameters, iMattox et al.l 119961 ) . We note that we require the predicted number of events 
from a source to be at least 10 over the full energy range, rejecting clusters of a few high- 
energy events without any low-energy counterpart. The same constraint was enforced for 
the 1FGL analysis. 

The analysis does have a number of important differences with respect to 1FGL: 



• The major change is that we switched from unbinned to binned likelihood (while still 
using gtlike or more precisely the pyLikelihood library in the Science Tools). The 
first reason for the change was to cap the computing time (which increases linearly 
with observing time in unbinned likelihood). The other important reason is that we 
discovered with simulations that the scale factors for the diffuse emission model terms 
returned by unbinned likelihood were significantly biased (overestimating the Galactic 
diffuse or isotropic diffuse intensity, whichever component was subdominant) whereas 
those returned by binned likelihood were not. In order to treat the Front and Back 
events in the analysis according to their separate PSFs we added the log£ computed 
separately for Front and Back events. The energy binning was set to 10 bins per 



5 See http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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decade. Rols are square for binned likelihood. We used the ARC projection with pixel 
size set to 0?1 for Front and 0?2 for Back events, in keeping with the high-energy PSF 
for each category. The sides of the Rols were defined by adding 7° on each side to 
the diameter of the central part where all source parameters are free. We note that 
the binned likelihood scheme is more conservative: in simulations comparable to the 
catalog depth (with or without sources) the significances of detections with unbinned 
likelihood tended to be around 1 a greater. This has important consequences for the 
number of sources in 2FGL (see § 14. 2p . 

• We took advantage of the fact that the localization procedure (§ I3.1.4P also provides 
a spectral fit to all sources. We used it as the starting point for the procedure using 
gtlike, rather than starting with all sources set to 0. 

• We did not use exactly the result of the previous iteration to start the next one, but 
applied a damping factor 5 (set to 0.1) to all parameters, defining the next starting 
point as P n +i = (1 _ 8)P n + SP n -\. It is a significant change because in all Rols the 
number of sources (outside the core of the Rol) which are considered but frozen is 
much larger than that of free sources. The damping procedure avoids overshooting 
and improves convergence. 

• Many bright sources are fitted with curved spectra instead of power-law. This is 
described in § 13.31 In addition to providing more detailed descriptions of those bright 
sources, it also improves the reliability of the procedure for neighboring sources. The 
reason is that it greatly reduces the spectral residuals, which otherwise might have 
been picked up by neighboring sources. That kind of transfer can be an issue at low 
energy where the PSF is very broad and cross-talk between sources in the likelihood 
analysis is strong. 

• We introduce the Earth limb component obtained in § 12.21 without any adjustment or 
free parameter in the likelihood analysis. 

Appendix [B] illustrates how well the full model (diffuse emission and individual sources) 
reproduces the 7-ray sky. 

3.3. Spectral Shapes 

The 1FGL catalog considered only power-law (PL) spectra. This was simple and ho- 
mogeneous, but not a good spectral representation of the bright sources, as could be easily 
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seen from comparing the power- law fits with the fluxes in bands (quantified by the Curva- 
tureJndex column in lAbdo et al.ll2010gl ). As the exposure accumulated, the discrepancies 
grew statistically larger, to the point where it could affect the global fit in an Rol, altering 
the spectra of neighboring sources in order to get a better overall spectral fit. For 2 years of 
data we had to allow for spectra that deviate from power laws. However increasing the num- 
ber of free parameters means finding the true best fit is more difficult, so we chose spectral 
shapes with only one additional free parameter. 

For the pulsars we chose exponentially cutoff power-l aws (hereafter PL ExpCutoff), 
which are a good representation of pulsar spectra in general flAbdo et al.ll2010tl ): 



dN 
dE 



cxp 



E-E c 



(1) 



This is just the product of power law and an exponential. The parameters are K, T (as 
in the power law) and the cutoff energy E c . E is a reference energy that we are free to 
choose for each source. The value of Eq started at 1 GeV but evolved separately for each 
source at each iteration as described below. All the known 7-ray pulsars with significant 
LAT pulsations were fitted with the PLExpCutoff representation. 

Other bright sources (mainly AGN ) are also not very well represented by power-law 
spectra. Analysis of the bright blazars (lAbdo et al.l l2010rl ) indicated that a broken power 
law was the best spectral representation. This however would add two free parameters 
and therefore was not stable enough for moderately bright sources. We adopted instead a 
log-normal representation (that we call LogParabola) which adds only one parameter while 
decreasing more smoothly at high energy than the PLExpCutoff form: 



dN 
dE 



K 



E 
En 



-a-piog{E/E ) 



(2) 



The parameters are K, a (spectral slope at E ) and the curvature (3, and E is an 
arbitrary reference energy that evolves for each source along the iterations. Negative 
(spectra curved upwards) were allowed, although we did not get any. 

In order to limit the number of free parameters, we did not fit every non-pulsar source as 
LogParabola, but only those in which the curvature was significant. In a procedure similar 
to that applied in the all-sky optimization for the source detection step (§ !3.1.2p . we assessed 
that significance for a given source by TS curve = 2 (log £ (LogParabola) — log £ (power-law)), 
where £ represents the likelihood function, changing only the spectral representation of 
that source and refitting all free parameters in the Rol. Since power-law is a special case 
of LogParabola (/3 = 0) and (3 = is inside the allowed interval we expect that TS curve is 
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distributed as x 2 with one degree of freedom. We switched to LogParabola if TS curve > 16, 
corresponding to 4 a significance for the curvature. All power-law sources were tested after 
each iteration, and we checked at the last iteration that TS curve for LogParabola sources 
was still > 16 (if it was not, the source was downgraded to power law and the Rol was 
refit). TScurve was computed for the LAT pulsars as well, but they were not downgraded to 
power-law if TS curve < 16. 

The extended sources (§ 13.4)) were handled on a case by case basis and fitted with either 
PLExpCutoff, LogParabola or power-law. 

The pivot energy E p (reported as Pivot_Energy) was computed as the energy at which 
the relative uncertainty on the differential flux K was minimal. This was done in the 
parabolic approximation using the covariance matrix between parameters. To improve the 
validity of the parabolic approximation, we changed the reference energy E used for fitting 
to E p after each iteration (with the same damping procedure as in § 13.21) . This ensured that 
at the end Eq was close enough to E p . The value of a (for LogParabola) depends on the 
reference energy, a(E p ) = a(E ) + 2(3\og(E p /E ). The uncertainties on K and a at E p were 
derived from the covariance matrix on the actual fitted parameters (relative to E ). The 
other parameters do not depend on the choice of E . 

In the catalog the differential flux K is reported as Flux_Density at the reference energy 
Eq = E p (where it is best determined). The low energy spectral index T (for PLExpCutoff) 
or the spectral slope Oi{E p ) (for LogParabola) are reported as SpectraLIndex. The cutoff 
energy E c is reported as Cutoff. The curvature (3 is reported as beta. For consistency with 
1FGL and in order to allow statistical comparisons between the power-law sources and the 
curved ones, we also report the spectral index of the best power-law fit as PowerLaw_Index 
for all sources. 

The fitted curvatures /3 sometimes tended to a large value, corresponding to very peaked 
spectra. There were cases (for example suspected millisecond pulsars) when this kind of 
spectrum could be real. However this occurred particularly in densely populated regions of 
the Galactic ridge, where the PSFs overlap and cross-talk between sources in the likelihood 
analysis is large at low energy. Even though one highly curved spectrum could lead to a 
better global fit for the Rol, it was not necessarily robust for that particular source, and in 
many cases we noted that the band fluxes (§ 13 . 5)) did not agree with the very curved fits. In 
order to avoid extreme cases, we enforced the condition < 1, corresponding to changing 
spectral slope by 2 log 10 = 4.6 over one decade. Whenever (3 reached 1 for a particular 
source, we fixed it to 1 and refitted in order to have a reasonable estimate of the errors 
on the other parameters. Sixty-four sources affected by this are flagged (see § 3.10 for a 
description of the flags used in the catalog). A similar difficulty occurred for 3 faint pulsars 
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in which the low energy index F tended to be very hard. We limited the values to T > 
and refitted with T fixed to when it was reached. Those 3 pulsars were flagged in the same 
way. Note that fixing one parameter tends to result in underestimating the errors on the 
photon and energy fluxes of those sources. 



3.4. Extended Sources 

In the analysis for the 1FGL catalog it became clear that a small number of sources were 
not properly modeled by a point source, leading to multiple detections being associated with 
the same source, e.g., the Large Magellanic Cloud (LMC). For the present analysis, twelve 
sources that have been shown to be extended in the LAT data were included as extended 
sources. The spatial templates were based on dedicated analysis of each source region, and 
have been normalized to contain the entire flux from the source (> 99% of the flux for 
unlimited spatial distributions such as 2-D Gaussians, 2DG). The spectral form chosen for 
each source is the closest of those used in the catalog analysis (see § 13.31) to the spectrum 
determined by the dedicated analysis^. 

The extended sources include seven supernova remnants (SNRs), two pulsar wind neb- 
ulae (PWNe), the LMC and the Small Magellanic Cloud (SMC), and the radio galaxy 
Centaurus A. Notes of interest for each source are provided below: 

• SMC - (2DG, PLExpCutoff) We modeled the SMC using a two-dimensional (2-D) 
Gaussian function with a width o = 0?9. While this is the best- fitting simple geom etric 



model, the morphology of the emission may be more complex (jAbdo et al.ll2010el ). 



LMC - (2x2DG, PLExpCutoff) This complex region, which accounted for five point 
sources in the 1FGL catalog, has been modeled as a c ombination of two 2 -D Gaussian 



profiles using the parameters specified in Table 3 of lAbdo et al.l (l2010ql ). The first, 
with a width of a = 1?2, represents emission from the entire galaxy. The second, 
with a width of o = 0?2, corresponds to the 7-ray bright region near 30 Doradus. 
Although this model provides a reasonable first order description of the 7-ray emission 
seen from the LMC, it is clear that this composite geo metric model is not sufficient to 
fully describe the complex morphology of the source JAbdo et al.ll2010qh . There are 



five sources in the 2FGL catalog that may be due to excess LMC emission after the 
fit, though two have blazar associations. 



6 The templates and spectral models will be made available through the Fermi Science Support Center. 
See Appendix ICl 
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IC 443 - (2DG, LogParabola) This SNR is modeled by a 2-D Gaussian profile with a 
width of a = 0?26. The log-parabola spectral fo rm most closely m atches the spectrum 
found for this source in the dedicated analysis (lAbdo et al.ll2010pl ). 



Vela X - (Disk, PL) We mode led Vela X using a simple disk with radius r = 0?88 
and a power law spectral form ( lAbdo et al.ll2010lf ). Since the Vela pulsar is spatially 
coincident with the Vela X PWN and significantly brighter, the detailed analysis was 
performed using the off-pulse events. For the catalog analysis it was necessary to fix 
the spectral parameters for the power law to the values determined by the off-pulse 
analysis. 

Centaurus A - (map, PL) This large radio galaxy has 7-ray emitting lobes that 
extend ~ 10° across the sky. The template used for this source originated from th e 
22 GHz WMAP image, and excludes a 1° region around the core (lAbdo et al.ll2010fl ). 
which is modeled separately as a point source in the catalog. The lobes are clearly 
resolved in the LAT. 

MSH 15—52 - (Disk, PL) This PWN is spatially coincident with the bright 7-ray 
pulsar PSR B1509— 58. The PWN was detected above 1 GeV, while the pulsar was 
detected only below 1 GeV by the LAT. We were able to investigate the PWN emission 
using events from all pulsar phases by excluding data below 1 GeV. T hat analysis 
showed that a uniform disk with radius r = 0?249 best fit the LAT data ( lAbdo et al. 



2010dl ). As with Vela X, the power-law spectral parameters for this source were fixed 
during the catalog analysis. 

W28 - (Disk, LogParabola) For W28, only the northern source at (R.A., Dec.) = 
(270?34, —23? 44) showed evidence for extension. We modeled this source using a disk 
with radius r = 0?39, the best-fit spatial model found by detailed analysis ( lAbdo et al. 



2010kJ). As with IC 443, a log-parabola spectral form fits the LAT data best. 



W30 - (Disk, LogParabola) The model for W30 uses a simple disk template centered 
at (R.A., Dec.) = (271?40, -21?63) with a radius r = 0?37. For the catalog analysis, 
a log-parabola spectral model best fits the source spectrum. 

HESS J1825-137 - (2DG, PL) This SNR is modeled with a 2-D Gaussian profile 
with a width of a = 0?56, which we found fit the source emission better than a disk. 
We tested a power-law spectrum both with and witho ut an exponential c utoff and 
found that the data was best fit by a simple power-law ( jGrondin et al.ll201ll ). 



W44 - (Ring, LogParabola) The template for the W44 SNR is an elliptical ring with 
axes (a, 6) inner = 0?22, 0?14, (a, 6) outer = 0?30, 0?19 and a position angle = 146° 
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counterclockwise from north (jAbdo et al.ll2010of ). Again, the best spectral model for 
the SNR is a log-parabola. 

W51C - (Disk, LogParabola) W51C is well r epresented by an el liptical disk with axes 
(a, b) = 0?40, 0?25 and a position angle 8 = 0° (jAbdo et al.ll2009bl ). using a log-parabola 
spectral form. 

Cygnus Loop - (Ring, PLExpCutoff) This relatively large SNR accounted for four 
sources in the 1FGL catalog. It is best represented by a ring located at (R.A., Dec.) = 
(312?75, 30?85) with an outer radius of r outer = 1?6 and an inner radius of r inner = 0?7 
jKatadri et alJboilh . 



Table [2] lists the source name, spatial template description, spectral form and the ref- 
erence for the dedicated analysis, where available. In the 2FGL catalog these sources are 
tabulated with the point sources, with the only distinction being that no position uncertain- 
ties are reported (see § 3.1.4). 



3.5. Flux Determination 

The source photon fluxes are reported in the 2FGL catalog in the same five energy 
bands (100 to 300 MeV; 300 MeV to 1 GeV; 1 to 3 GeV; 3 to 10 GeV; 10 to 100 GeV) 
as in 1FGL. The fluxes were obtained by freezing the spectral index to that obtained in 
the fit over the full range and adjusting the normalization in each spectral band. For the 
curved spectra (§ 13. 3j) the spectral index in a band was set to the local spectral slope at the 
logarithmic mid-point of the band y/E n E n+1 , restricted to be in the interval [0,5]. We used 
binned likelihood in all bands, but contrary to § 13.21 we did not distinguish Front and Back 
events. The pixel sizes in each band were 0?3, 0?2, 0?15, 0?1, 0?1 decreasing in size with 
energy as the PSF improves. 

The procedure for reporting either a measurement or an upper limit is the same as for 
the 1FGL catalog. For bands where the source was too weak to be detected, those with Test 
Statistic in the band TSi < 10 or relative uncertainty on the flux AFj/Fj > 0.5, 2 a upper 
limits we re calculated, FY L . Two methods were used, the profile and Bayesian methods. In 



the first (IRolke et al.ll2005l ). which is used when 1 < TS < 10, the profile likelihood function 
log£(Fj), is assumed to be distributed as x 2 /2 and the upper limit corresponds to th e point 



wher e log£(Fj) decreases by 2 from its maximum value. In the Bayesian method (IHelene 



19831 ). which is used when TS < 1, the limit is found by integrating C(Fi) from up to the 



flux that encompasses 95% of the posterior probability. With the probability chosen in this 
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Fig. 4. — A representative map of the TS residuals along the Galactic ridge, corresponding 
to the final iteration. There are several clusters that could have generated seeds, but did not 
appear to be isolated point sources under visual inspection. 



Table 2. Extended sources used in the 2FGL analysis 



2FGL Name 


Extended Source 


Spatial Form 


Spectral Form 


Reference 


2FGL 


J0059.0-7242c 


SMC 


2D Gaussian 


Exp Cutoff PL 


Abdo et al. (2010e1 


2FGL 


J0526.6-6825e 


LMC 


2D Gaussian a 


Exp Cutoff PL 


Abdo et al. (2010a 1 


2FGL 


J0617.2+2234c 


IC 443 


2D Gaussian 


Log Parabola 


Abdo et al. (2010d1 


2FGL 


J0833.1-4511e 


Vela X 


Disk 


Power Law 


Abdo et al. ( 20101) 


2FGL 


J1324.0-4330e 


Centaurus A (lobes) 


Contour Map 


Power Law 


Abdo et al. (2010f) 


2FGL 


J1514.0-5915e 


MSH 15-52 


Disk 


Power Law 


Abdo et al. f2010d) 


2FGL 


J1801.3-2326c 


W28 


Disk 


Log Parabola 


Abdo et al. (2010k1 


2FGL 


J1805.6-2136c 


W30 


Disk 


Log Parabola 




2FGL 


J1824.5- 1351c 


HESS J1825-137 


2D Gaussian 


Power Law 


Grondin et al. (20111 


2FGL 


J1855.9+0121c 


W44 


Ring 


Log Parabola 


Abdo et al. C2010o1 


2FGL 


J1923.2+1408e 


W51C 


Disk 


Log Parabola 


Abdo et al. (2009b1 


2FGL 


J2051.0+3040e 


Cygnus Loop 


Ring 


Exp Cutoff PL 


Kataairi et al. ("20111 



a To fit the LMC we used a combination of two 2D Gaussian spatial templates. 



Note. — Twelve 2FGL sources that have been modeled as extended sources. More detail regarding the 
parameters used in the analysis can be found in the text. The publications describing the detailed analysis 
for W30 is still in preparation. 



-27- 



way the upper limits calculated with each method are similar for sources with TS = 1. The 
2 a upper limit is then reported in the flux column and the uncertainty is set to 0. 

In the 1FGL catalog the photon flux betwee n 1 and 100 GeV an d the energy flux between 



100 MeV and 100 GeV (F 35 and S 25 in Table 151 lAbdo et al.ll2010gf > were estimated from the 
sum of band fluxes because the result of the fit over the full band was biased by the power-law 
approximation and was inconsistent with the sum of band fluxes for the bright sources. In 
the 2FGL catalog analysis the curved spectral shapes are precise enough to overcome that 
limitation (Fig. [SJ. The main advantage of the full spectral fit is that it is statistically more 
precise because it incorporates the (reasonable) constraint that the spectral shape should be 
smoothly varying with energy. Even using the newer data set (with larger effective area at 
low energy), the relative uncertainties in the lower energy bands tend to be very large. The 
relative uncertainty on the full photon flux between 100 MeV and 100 GeV (-F25, dominated 
by low energy) is much larger than that on F35 or S25 (23% vs 15% and 14% respectively for 
a TS = 100 source with spectral index 2.2) and strongly depends on spectral index (whereas 
that on F35 does not). So we do not report the photon flux over the full band in 2FGL. 
We report F 35 and S25, as in 1FGL, but estimated from the fit over the full band. For 
comparison, the relative uncertainties on estimates of F 35 and S25 from the sum of bands 
(as in 1FGL) are 20% for the same typical source. The procedure for reporting upper limits 
described above applies to F 35 and S25 as well. Five sources (4 very hard and 1 very soft) 
have relative uncertainty on F35 larger than 0.5. The faintest of those 5 also has relative 
uncertainty on S25 larger than 0.5. 

We show the photon and energy flux distributions for the 2FGL sources in two different 
ways in Figures [6] and [7J Figure UJ shows that t he range of energy fl uxes among the 2FGL 



sources is greater than 3 decades. Figure 20 of lAbdo et al.l (l2010gl ) was the same plot as 
Figure but on the photon flux between 100 MeV and 100 GeV. The detection threshold on 
the photon flux over the full band depends sensitively on the spectral index of the source. 
Building a flux-limited sample on that quantity required raising the minimum flux to the 
detection threshold for soft sources and resulted in discarding most of the hard sources. 
The photon flux above 1 GeV (or the energy flux), which we show in these figures, is more 
appropriate to build a flux-limited sample because it discards few sources. 

Figures [BJ [9j and [10] show examples of the band fluxes, with the best fit over the full 
range overlaid. From this kind of plot one may build a spectral fit quality indicator similar 
to the Curvature_Index of 1FGL. 

C V (fl-*f) 2 (3) 
where i runs over all bands and Ff l is the flux predicted in that band from the spectral fit to 
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Fig. 5. — Comparison of estimates of the energy flux from 100 MeV to 100 GeV S25 from 
the sum of bands (abscissa) and the fit to the full band (ordinate). No obvious bias can be 
observed. 
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Fig. 6. — Distribution of sources in 2FGL excluding the Galactic plane in the spectral index 
- photon flux plane. The spectral index is the effective PowerLaw_Index (power-law fit even 
for curved sources). The photon flux is between 1 and 100 GeV {F 35 ). The low flux threshold 
is quite sharp arou nd 4 x 10~ 1 ph cm " 2 s _1 . The full line shows the expected threshold 
following App. A of lAbdo et al.l (]2010gM accounting for the average confusion, and the dashed 
line for an isolated source. 



-30 - 



1000.0 FT 



100.0 - 



10.0 - 



CD 
42 



3 

z 



1.0 - 



0.1 I I , , , , ! I , I I I I ! 

jq-12 jQ-11 jQ-10 1Q -9 jQ-l 

Knergy flux in erg/cm2/s 



Fig. 7. — Distribution of all sources in 2FGL with respect to log(Energy flux). The low flux 
threshold is quite sharp around 5 x 10~ 12 erg cm -2 s _1 , indicating that the TS cut that 
is applied is not too far from a cut on the energy flux S 2 5 over the full band (100 MeV to 
100 GeV). 
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Fig. 8. — Spectrum of a faint AGN, as an example of a power-law spectrum. The fit over 
the full band (dashed line) is overlaid over the five band fluxes converted to vF v units. The 
grey shaded area (butterfly) shows the formal 1 a statistical error on log (differential flux) as 
a function of energy, obtained using the covariance matrix involving the parameters of that 
particular source. The upper limits (here the lowest-energy and highest-energy bands) are 2 
a. 
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Fig. 9. — Spectrum of the pulsar in CTAl, as an example of an exponentially cutoff spectrum. 
See Figure M for details. 
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Fig. 10. — Spectrum of the bright AGN 4C +21.35, as an example of a LogParabola spec- 
trum. See Figure M for details. 
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the full band, ft reflects the systematic uncertainty on effective area (§ 13. Tj) . They were set 
to O.f , 0.05, 0.05, 0.08, 0.1 in our five bands. Since, in 2FGL, curvature is accounted for in 
the spectral shape, the interpretation of that quantity is now whether the proposed spectral 
shape agrees well with the band fluxes or not. We did not report that in the table, but we 
set a flag (Flag 10 of Table EJ) whenever C sys t > 16.3, corresponding to a probability of 10~ 3 
assuming a \ 2 distribution with 3 degrees of freedom (5 — 2, since the majority of sources 
are fitted with power-law spectra which have 2 free parameters). Thirty-three sources are 
flagged in this way, including the two brightest pulsars (Geminga and Vela) whose spectrum 
does not decrease as fast as a simple PLExpCutoff. 

A few percent error in the effective area calibration as a function of energy may result in 
an incorrect report of significant curvature for very bright sources. There is no obvious rig- 
orous way to enter systematic uncertainties in the TS curve calculation (§ 13.31) . In order to do 
that approximately, we note that TS curve is an improved estimator of how much the spectrum 
deviates from a power-law. The analog of TS curve in 1FGL was C^ syst , applying Eq. [3] to the 
power-law fit with no // el term (TS curve is a purely statistical quantity). We can compare 
Cnrayst w ith the same quantity Cf^ t obtained with the // el term (Curvature_Index of 1FGL). 
Their ratio is a measure of how much the systematic uncertainties reduced Curvature_Index. 
We can then apply that same ratio to TS curve and we report in the catalog Signif _Curve = 



We consider that sources with Signif _Curve > 4 are significantly curved. The conse- 
quence of introducing the systematic uncertainties is that 40 sources in the catalog have a 
LogParabola spectrum because TS curve > 16 (§ 13. 3j) even though Signif _Curve < 4. We 
do not claim that the curvature is real for those sources, even though it is statistically 
significant. 



Temporal variability is relatively common in 7-ray sources and provides a powerful 
tool to associate them definitively with objects known at other wavelengths and to study 
the physical processes powering them. We present a light curve for each source in the 
catalog, produced by dividing the data into approximately monthly time bins and applying 
the likelihood analysis procedure to each. The details of the light curve analysis and how 
the results are presented are summarized below: 

• There are 24 time bins, starting at the beginning of the data set, approximately 



54682.66 MJD (§ 12. ip . The first 23 bins have durations of 30.37 days; the final has 




converting to a units. 



3.6. Variability 
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a duration of 27.88 days. The first 11 time bins correspond exactly to those of 1FGL. 

• The parameters describing the spectral shapes of the sources in the Rol are fixed in the 
light curve calculation. Only the normalizations of the source of interest, the diffuse 
backgrounds, and bright and nearby catalog sources (see section 13. 2j) are allowed to 
vary. We use binned likelihood, but do not distinguish Front and Back events. The 
pixel size is set to 0?2. 

• For each time bin, the photon flux over the full energy range (100 MeV to 100 GeV), 
Fi, its error, AFj and the detection significance, TSi, are presented in the catalog. 
With the spectral shape of each source frozen in the light curve analysis, the relative 
uncertainty on F 25 is the same as that of F 35 and S25, and it is reasonable to present 
the photon flux over the full energy range in this case. 

• For time bins where the source is too weak to be detected, those with TSi < 10 or 
AFi/Fi > 0.5, 95% upper limits Ff L are calculated following the same method as in 
§ 13.51 A fraction of those have flux exactly equal to 0, because the Poisson likelihood 
framework that we use does not accept negative flux values. 

• In the case of an upper limit, the best-fit flux value is given in the catalog, and the 
error is replaced by 0.5(F i ' 7i — Fi). This allows bands with upper limits to be treated 
consistently with the other bands while preserving enough information to extract the 
upper limits. The FITS version of the catalog^ has a flag column to indicate when an 
entry in a flux history is an upper limit. Please note that this is a different convention 
to that used to report flux upper limits for the energy bands (§ 3.5). See Appendix C 
for more information. 

• A total of 340 sources have only upper limits on monthly timescales. These sources 
have an average integrated significance over the full 2-year data set of 5.3 a. At the 
opposite extreme, 94 sources are detected significantly in every one of the time periods. 

To test for variability in each source we construct a variability index from the value of 
the likelihood in the null hypothesis, that the source flux is constant across the full 2-year 
period, and the value under the alternate hypothesis where the flux in each bin is optimized: 

TS var = 2 [log C({F t }) - hgC(Fconst)} = 2 ^ [log d{Fi) - log A(F Const )] = 2 ^ V? (4) 

i i 



7 The FITS version of the catalog is available through the Fermi Science Support Center. See Appendix [Cl 
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where the log likelihood for the full time period, log C({F{}), can be expressed as a sum 
of terms for the individual time bands, log£j. If the null hypothesis is correct TS var is 
distributed as x 2 with 23 degrees of freedom, and a value of TS var > 41.6 is used to identify 
variable sources at a 99% confidence level. For most sources the value for F Con st is close 
to the value derived from the likelihood analysis of the full time period, although strong 
variability in nearby background sources can cause to them to differ in some cases. The light 
curve for PKS 1510—089, a bright blazar, is shown in Figure [TT1 This source is easily flagged 
as variable, with TS var = 6406. 

Upper limits calculated through the profile method are handled naturally in the variabil- 
ity index procedure described above, but those calculated using the Bayesian method would 
have to included in an ad hoc manner. Instead, when calculating the variability index, the 
results of the profile method are used for all upper limits. 

As in 1FGL, the brightest pulsars detected by the LAT are flagged as being variable with 
this procedure. This apparent variability is caused by systematic errors in the calculation 
of the source exposure, resulting from small inaccuracies in the dependence of the IRFs on 
the source viewing angle, coupled with changes in the observing profile as the orbit of the 
spacecraft precesses. We introduce a correction factor to account for these errors, and fix 
the size of this correction such that the bright pulsars are steady. Specifically, we scale each 
V? in the summation of TS var by a factor which combines the error on the flux each time 
bin in quadrature with a fixed fraction of the overall flux, 

A P 2 

A value of / = 0.02, i.e. a 2% systematic correction factor, was found sufficient such that only 
PSR J1741— 2054 remains (marginally) above threshold among the LAT pulsars, excluding 



Tavani et al. 


2011; 


Abdo et al. 


2011c) 



required in 1FGL, the improvement resulting from the higher-fidelity IRFs used in this work. 
This systematic error component is included in the flux errors reported in the catalog FITS 
file. Figure [T2l shows the light curve for the pulsar Geminga (PSR J0633+1746), one of the 
brightest non-variable sources in 2FGL. 

The Sun is a bright, extended source of 7 rays, both from cosmic-ray interactions in 
its outer atmosphere and from IC scattering of cosmic-ray electr ons on the solar rad iation 



field, which produces an extended 7-ray halo around the Sun (lAbdo et al.l l2011bl ). For 
sources close to the ecliptic, solar conjunctions can lead to significant enhancements of the 
flux detected during the time periods when the Sun is closer than approximately 2?5 to the 
source. Sources for which a large fraction of the total detection significance comes during 
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Fig. 11. — Light curve for the bright blazar PKS 1510—089 in the full energy range (100 MeV 
to 100 GeV). The dashed line depicts the average flux from the analysis of the full 24- month 
dataset. 



-38 - 



2FGL J0633.9+1 746 - PSR J0633+1 746, Geminga 



4.45 




3.95 



54700 54800 54900 55000 55100 55200 55300 55400 

MJD [days] 



Fig. 12. — Light curve for the bright pulsar Geminga (100 MeV to 100 GeV). The gray band 
depicts the size of the 2% systematic correction applied to the calculation of the variability 
index. The error bars on the flux points show the statistical errors only. 



-39 - 



such periods are nagged as suspicious in the catalog. The light curve for such a source, 
2FG L J2124.0— 1513, is shown in Figure [13j The Moon is comparably bright to the Sun in 7 
rays (lAbdo et al.ll2012al ) and lunar conjunctions also potentially affect the fluxes measured 
from LAT sources, but the higher apparent speed of the Moon and parallax from the motion 
of the spacecraft mean that such close conjunctions are brief, and the precession of the orbit 
of the Moon means that for any given source conjunctions are less frequent than the ~28 d 
period of the orbit. In a ddition, the 7-ray e mission of the Moon does not include a bright, 
extended IC component ( lAbdo et al.ll2012ai ). Hence, we do not attempt to identify sources 
which may be affected by the Moon nor to flag time periods where lunar conjunctions occur. 



Light curves for all 2FGL sources are available from the Fermi Science Support Center. 



3.7. Limitations and Systematic Uncertainties 

A limitation for the catalog analysis is source confusion. (The related issue of systematics 
for localization is discussed in § 13.1.41 ) Confusion is of course strong in the inner Galaxy, 
where the source density is very high, but it is also a significant issue elsewhere. The average 
distance between sources outside the Galactic plane is 2?8 (it was 3° in 1FGL), to be compared 
with a per photon containment radius r^g = 0?8 at 1 GeV where the sensitivity is best. The 
ratio between these numbe rs is not large eno ugh that confusion can be neglected. As for 



the 1FGL catalog analysis (lAbdo et al.ll2010gl ) we study source confusion by evaluating the 



distribution of distances between each source and its nearest neighbor (D n ) in the area of 
the sky where the source density is approximately uniform, i.e., outside t he Gal actic plane. 



This is shown in Figure [HJ to be compared with Figure 9 of lAbdo et al.l ( l2010gt ) which also 
details the expected distribution. The histogram still falls off toward D n = 0, but follows 
the expected distribution down to 1° or so instead of 1?5 in 1FGL. We estimate that some 
43 sources within 1° of another one were missed because of confusion (to be compared with 
the 1319 sources observed at |6| > 10°). This means that the fraction of missed sources 
decreased from 7.7% in the 1FGL analysis to 3.3% for 2FGL. This attests to the progress 
made in the detection process (§ 13. ip . 

An important issue for the evaluation of spectra is the systematic uncertainties of the 
effective area of the instrument. Compared to the 1FGL instrument response functions 
(P6_V3), the current P7_V6 response functions have somewhat reduced systematic uncer- 
tainties. The current estimate of the remaining systematic uncertainty is 10% at 100 MeV, 



decrea sing to 5% at 560 MeV and increasing to 10% at 10 GeV and above (lAbdo et al. 



2012bl ). This uncertainty applies uniformly to all sources. Our relative errors (comparing 



one source to another or the same source as a function of time) are much smaller, as indicated 
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Fig. 13.— Light curve for the unassociated source 2FGL J2124. 0-1513 (100 MeV to 
100 GeV). Time periods in which the sun is closer than 2?5 to the source are marked with yel- 
low vertical bands. In this case, a large fraction of the detection significance is accumulated 
during these periods, and the source is flagged as suspicious in the catalog. 
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Fig. 14. — Distribution of the distances D n to the nearest neighbors of all detected sources 
at \b\ > 10°. The number of entries is divided by 27iD n AD n in which AD n is the distance 
bin, in order to eliminate the 2-dimensional geometry. The overlaid curve is the expected 
Gaussian distribution for a uniform distribution of sources with no confusion. 
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in §EU 

The model of diffuse emission is the other important source of uncertainties. Contrary 
to the effective area, it does not affect all sources equally: its effects are smaller outside the 
Galactic plane where the diffuse emission is fainter and varying on larger angular scales. It 
is also less of a concern in the high energy bands (> 3 GeV) where the core of the PSF is 
narrow enough that the sources dominate the background under the PSF. But it is a serious 
concern inside the Galactic plane in the low-energy bands (< I GeV) and particularly inside 
the Galactic ridge (|/| < 60°) where the diffuse emission is strongest and very structured, 
following the molecular cloud distribution. It is not easy to assess precisely how large the 
uncertainties are, because they relate to uncertainties in the distributions of interstellar gas, 
the interstellar radiation field, and cosmic rays, which depend in detail on position on the 
sky. We discuss the Galactic ridge more specifically in § 13.91 

For an automatic assessment we have tried re-extracting the source locations and fluxes 
assuming the same diffuse model that we used for fFGL, and also the same event selection 
as in 1FGL but with improved calibration (P6_vff). The results show that the systematic 
uncertainty more or less follows the statistical one (i.e., it is larger for fainter sources in 
relative terms) and is of the same order. More precisely, the dispersions of flux and spectral 
index are 0.8 o at |6| > 10°, and f.3 a at |6| < 10°. We have not increased the errors 
accordingly because this older model does not fit the data as well as the newer one. From 
that point of view we may expect this estimate to be an upper limit. On the other hand 
both models rely on nearly the same set of H I and CO maps of the gas in the interstellar 
medium, which we know are an imperfect representation of the mass. That is, potentially 
large systematic uncertainties are not accounted for by the comparison. So we present the 
figures as qualitative estimates. We also use the same comparison to flag outliers as suspect 

(§ED5. 

Finally, we note that handling Front and Back events separately for the significance 
and spectral shape computation (§ 13. 2j) introduces another approximation. Because the free 
parameters are the same for both categories of events, this amounts to assuming that the 
isotropic diffuse emission is the same for Front and Back events. This is actually not true 
because it contains internal background that is larger for Back events (see §I2])- This effect 
is significant only below 400 MeV (§ 12. 2p . and so the consequence is an underestimate of the 
low-energy flux, which results in a systematic hardening of the measured power-law spectral 
index but which is nearly always less than its statistical uncertainty. Thus in terms of the 
absolute change in spectral index, the effects are greatest for soft sources. 



3.8. Point Sources and Extended Sources 



Except for the diffuse emission and the 12 sources explicitly considered as spatially 
extended, all sources in the catalog are assumed to be point-like. Just as the modeling of 
the diffuse emission can affect the properties of point sources (as discussed in the previous 
section), the treatment of known or unknown extended sources can similarly influence the 
analysis of nearby point sources. This influence can be felt in three ways: 



1. The modeling of an extended source is limited by the detailed knowledge of the 7-ray 
emissivity of the source as a function of position on the sky. As noted in § 3.4, the 
modeling for the catalog was done using largely geometric functions. The true distribu- 
tion can have residual excesses that the catalog analysis then treats as point sources. 
Examples are the sources near the Large Magellanic Cloud: 2FGL J0451.8— 7011, 
2FGL J0455.8-6920, 2FGL J0532.5-7223, 2FGL J0533.3-6651, 2FGL J0601. 1-7037 
Although some of these may be unrelated to the LMC itself (two have blazar associa- 
tions), some may be residuals from the modeling. Sources close to any of the extended 
sources should be treated warily in detailed analysis of such regions. 

2. Some known or likely extended sources are not among the 12 that were modeled for the 
catalog analysis, having been recognized and measured only after the catalog analysis 
was largely complete. In such cases, the catalog analysis finds one or more point 
sources at or near the possible extended source. Two clear examples are supernova 
remnants. RX J1713.7-3946 is represented in the catalog by 2FG L J1712.4-3941, bu t 



recent analysis has shown this SNR to be an extended GeV source ( lAbdo et al.ll2011dl ). 
In Table EH RX J0852.0-4622 shows four associated 2FGL sources: J0848.5-4535, 
J0851. 7-4635, J0853.5-4711, and J08 55.4-4625. All of these are likely part of the 



spatially extended supernova remnant (ITanaka et al. 2011 ). Other clusters of sources 



in Table [TT] may indicate yet-unresolved extended objects. As longer exposures with 
the LAT collect more of the highest-energy photons with the best angular resolution, 
additional spatial structure will be revealed in the data. 

3. A spectral bias can be introduced if an extended source is analyzed as if it were a 
point source. In such cases the calculated spectrum is likely to be softer than the true 
spectrum. At higher energies where the LAT PSF is closer to the size of the extended 
source, the extension will cause such photons to be lost. 
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3.9. Sources Toward Local Interstellar Clouds and the Galactic Ridge 

The interstellar part of the model for diffuse emission of the Galaxy has greatly improved 
since the 1FGL catalog analysis, in particular in angular resolution (§ 12 . 2 j) . However, the 
use of large-scale rings in the Milky Way and of a single ring in the solar neighborhood 
(containing most of the gas-related diffuse emission off the Galactic plane) does not allow 
for small-scale variations in the gas and dust properties used to derive the target mass for 
cosmic rays, or in the cosmic-ray spectrum itself. The optical depth correction applied in 
deriving column densities from the H I 21-cm line observations also is necessarily based on 
the approximation that the gas has a uniform spin temperature across the Galaxy. As a 
result, extended and structured excesses of 7 radiation are present above the diffuse model. 
The renormalization of the diffuse model within each Rol lessens, but cannot always remove, 
the impact of the diffuse excesses. Point sources detected in the regions of the excesses can 
be formally very significant but are necessarily suspect (see Fig. ITS"]) . 



One particularly difficult region is the local arm tangent in Cygnus. The region 75° < 
/ < 85°, |fe| < 15° contains 18 2FGL sources, 13 without firm identifications. Of these 13, 
we find that 8 sources, with detection significances ranging up to 13 a have properties that 
are especially sensitive to the diffuse emission model (see § I3.10[ in particular Flags 1-4 in 
Table El) . The diffuse em ission in Cygnus has recently been the subject of a detailed study 



( jAckermann et al.l l2012bl ) and evidence suggests that some of the exces s diffuse emission 



is due to an extended cocoon of unusually hard-spectrum cosmic rays ( jAckermann et al. 



2011al ). 



We have used a dust reddening map to trace substantial amounts of dark gas in addition 



to the atomic an d molecular gas seen in H I and CO emission lines ( IGrenier et al.l 12005 



Ade et al.ll2011b[ ). This made essential improvements over wide regions from low to medium 



l atitudes, but inaccu racies in the infrared color corrections used to build the reddening map 



( iSchlegel et al.lll998l ) can cause diffuse residuals toward bright H II regions or stellar clusters 



by artifici ally lowering the g as column densities measured in their directions (see Figs. 10 



and 11 of lAbdo et al.ll2010gl ). There are fewer such artifacts in 2FGL than in 1FGL, but 
examples can be found in Orion, Taurus, and near the source LSI +61 303; see also Figure 
[TBI Another known limitation of the diffuse model relates to the optical thickness of the CO 
lines and the saturation of the CO intensity toward very dense clouds. Since stellar clusters 
are born in the clouds, both CO saturation and dust temperature corrections may cooperate 
to under-predict the gas mass in dense molecular clouds. Self-absorption of the H I lines also 
leads to under-predicted column densities in the dense atomic phase. These limitations are 
particularly relevant at low latitudes, in the inner Galaxy or toward the tangent directions 
of the Galactic spiral arms. 
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Fig. 15. — From top to bottom: the CO contribution to the interstellar photon counts, 
the total interstellar photon counts, and the photon residual counts above the model for 
diffuse 7-ray emission, all in the 1-11 GeV energy band. The circles mark the effective 50% 
containment radii of the 2FGL sources for the 1-10 GeV band. The 'c' sources are crossed. 
The square notes an identified source. The photon residual map has been smoothed for 
display with a o = 0?125 Gaussian. The 2FGLc sources seen above the Galactic plane, with 
TS ranging from 26 to 75, follow an extended and clumpy excess of interstellar emission 

We have inspected all the 2FGL sources to search for potential problems with the 
underlying diffuse model. It is unlikely that sources with very high TS can be diffuse 
excesses. Based on the examination of the sources toward Cygnus, Orion and other nearby 
clouds, as well as the 1FGL sources with the 'c' designation that are not confirmed in 2FGL, 
we tentatively consider that sources with TS > 200, 130, or 80 are unlikely to be diffuse 
features depending on the intensity of the diffuse background (respectively when the photon 
count per pixel N bkgd , integrated from 589 MeV to 11.4 GeV in the diffuse model cube, 
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Fig. 16. — From top to bottom: the absolute value of the dust negative residual photon 
counts incorporated in the diffuse model, the total interstellar photon counts, and the photon 
residual counts above the diffuse model, all in the 1-11 GeV energy band. The circles mark 
the effective 50% containment radii of the 2FGL sources for the 1-10 GeV band. The 'c' 
sources are crossed. The photon residual map has been smoothed for display with a a = 0?125 
Gaussian. The 2FGL 'c' sources are distributed along the rim of a large H II region where 
the dust temperature correction led to an overestimate of the dust column densities in the 
ionized gas. The negative dust residuals have artificially reduced the diffuse 7-ray intensity 
in these directions. 



without the isotropic contribution, is N bkgd > 100, 60 < N bkgd < 100, or N bkgd < 60). 

Given the large change in the width of the PSF across the LAT energy band, we com- 
puted the effective 50% containment radius for each source from its best-fit spectrum. We 
overlaid these on predicted photon count maps from the Galactic diffuse model, both for the 
total emission and for the individual gas components in each phase, in seven energy bands 
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(the five energy bands of the catalog, plus the integral 0.5-10 GeV and 1-10 GeV bands). We 
also compared photon residual maps (data minus model) in the same energy bands against 
the predicted counts maps for the individual gas components. We also took into account the 
TS values reached in the five catalog energy bands and the spectral index of each source. Off 
the Galactic plane, we flagged (Flag 6 of Table |3]) unassociated sources coinciding with dust 
temperature or dense CO defects, or concurrent with extended residuals that followed inter- 
stellar features (as in Figure ITS"]) . Sources with TS larger than the background-dependent 
threshold quoted above or with a spectral index T < 2 were not flagged. In the Galactic 
plane (i.e. at |6| < 2° for \l\ < 70°, or |6| < 3° at higher longitudes), we flagged two types of 
sources: (i) unassociated sources with overlapping 50% containment radii above 500 MeV, 
unless their TS exceeded the background-dependent threshold or their spectral index were 
< 2; (ii) low-significance unassociated sources with TS < 80 for N bkgd > 160, unless their 
spectral indices were < 2. This strategy ensured that most of the Galactic ridge sources, 
which are closely packed together to make up for the extended photon residuals along the 
plane, are flagged, but it leaves all the identified and associated, intense, and hard sources 
out of the systematic ridge flag we had used in 1FGL. Every source was then manually 
checked with the same set of maps as for the work at higher latitude. 

We have added the designator 'c' to the names of the flagged sources to indicate that 
they are to be considered as potentially confused with interstellar emission. Their position, 
emission characteristics, or even existence may not be reliable. The 'c' designator applies to 
162 sources in the 2FGL catalog. 

3.10. Analysis Flags 

As in 1FGL we identified a number of conditions that can shed doubt on a source. They 
are described in Table [3] In the FITS version of the catalog, these flags are summarized in 
a single integer column (Flags). Each condition is indicated by one bit among the 16 bits 
forming Flags. The bit is raised (set to 1) in the dubious case, so that sources without any 
warning sign have Flags = 0. 

Flags 1 to 9 have similar intent as in 1FGL, but differ in detail: 

• In Flag 4, we reduced the threshold on source to background ratio to 20%, because the 
diffuse model has improved. 

• The distances triggering Flag 5 have changed because the PSF knowledge has improved. 
The core of the PSF at low energy is actually better than the P6V3 estimate used in 
1FGL, so the critical distance is lower at low energy. On the other hand the measured 
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Table 3. Definitions of the Analysis Flags 



Flag a Meaning 

1 Source with TS > 35 which went to TS < 25 when changing the diffuse model 
(§ I3.7|) . Note that sources with TS < 35 are not flagged with this bit because 
normal statistical fluctuations can push them to TS < 25. 

2 Moved beyond its 95% error ellipse when changing the diffuse model. 

3 Flux (> 1 GeV) or energy flux (> 100 MeV) changed by more than 3 a when 
changing the diffuse model. Requires also that the flux change by more than 
35% (to not flag strong sources). 

4 Source-to-background ratio less than 20% in highest band in which TS > 25. 
Background is integrated over 7rr| 8 or 1 square degree, whichever is smaller. 

5 Closer than # re f from a brighter neighbor. # re f is defined in highest band in 
which source TS > 25, or the band with highest TS if all are < 25. 9 rc{ is set 
to 2? 17 (FWHM) below 300 MeV, 1?38 between 300 MeV and 1 GeV, 0?87 
between 1 GeV and 3 GeV, 0?67 between 3 and 10 GeV and 0?45 above 

10 GeV (2r 68 ). 

6 On top of an interstellar gas clump or small-scale defect in the model of 
diffuse emission; equivalent to the 'c' designator in the source name (§ 13.91) . 

7 Not used. 

8 Inconsistent position determination (§ 13. 1.4[) : best position from optimization 
outside the 1 o (39% in 2D) contour from the TS map. 

9 Elliptical quality > 4 in 'pointlike (i.e., TS contour does not look elliptical). 

10 Spectral Fit Quality > 16.3 (EqJBJ. 

11 Possibly due to the Sun (§ 13.61) . 

12 Highly curved spectrum; LogParabola ft fixed to 1 or PLExpCutoff 
SpectraLIndex fixed to (see § I3.3j) . 



a In the FITS version the values are encoded as individual bits in a single column, with 
Flag n having value 2 <n_1 < ) . For information about the FITS version of the table see Table [T21 
in Appendix O 



in-fligh t PSF at high energy is much broader than the P6V3 estimate (lAbdo et al. 
2009hl ). so the critical distance is about twice as great than for the 1FGL analysis 
above 10 GeV. 



• We do not use gtfindsrc in 2FGL because it is based on unbinned likelihood. Therefore 
Flag 7 is not used. 

• Flag 8 compares the best position obtained from direct optimization with the contours 
extracted from the TS maps. 

• The threshold for Flag 9 on elliptical quality was decreased to 4. The improved local- 
ization procedure allowed being a little more stringent here. 

Flags 10, 11, and 12 are new. Figure [17] shows the distribution on the sky of flagged 2FGL 
sources. 

4. The 2FGL Catalog 

The basic description of the 2FGL catalog is in § 14.11 including a listing of the main 
table contents and some of the primary properties of the sources in the catalog. We present 
a detailed comparison of the 2FGL catalog with the 1FGL catalog in § 14.21 

4.1. Catalog Description 

Table S] is the catalog, with information for each of the 1873 sources; see Table [5] for 
descriptions of the columns. The source designation is 2FGL JHHMM . m+DDMM where the 2 
indicates that this is the second LAT catalog, FGL represents Fermi Gamma-ray LAT. Sources 
close to the Galactic ridge and some nearby interstellar cloud complexes are assigned names 
of the form 2FGL JHHMM. m+DDMMc, where the c indicates that caution should be used in 
interpreting or analyzing these sources. Errors in the model of interstellar diffuse emission, 
or an unusually high density of sources, are likely to affect the measured properties or even 
existence of these 162 sources (see § 13. 9j) . In addition a set of analysis flags has been defined 
to indicate sources with unusual or potentially problematic characteristics (see § 13 .91) . The 
'c' designator is encoded as one of these flags. An additional 315 sources have one or more of 
the other analysis flags set. The 12 sources that were modeled as extended for 2FGL (§ I3.4[) 
are singled out by an e at the ends of their names. 
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Fig. 17. — Full sky map (top) and blow-up of the inner Galactic region (bottom) showing 
flagged sources by source class. Sources potentially confused with diffuse emission, i.e., those 
with a 'c' designator in their names (and for which Flag 6 is set) are shown in red; those 
with any other flag set are shown in blue. Sources with no flag set are shown as small dots. 



-51 - 



The designations of the classes that we use to categorize the 2FGL sources are listed in 
Table |6] along with the numbers of sources assigned to each class. We distinguish between 
associated and identified sources, with associations depending primarily on close positional 
correspondence (see § 15. 2p and identifications requiring measurement of correlated variability 
at other wavelengths or characterization of the 2FGL source by its angular extent (see § 15. ip . 
In the cases of multiple associations with a 2FGL source, we adopt the single association 
that is statistically most likely to be true if it is above the confidence threshold (see § 15. 2j) ; 
the one exception is the Crab pulsar and PWN, which are listed as being associated with the 
same 2FGL source (see 15.11) . Sources associated with SNRs are often also associated with 
PWNs and pulsars, and the SNRs themselves are often not point-like. We do not attempt 
to distinguish among the possible classifications and instead list plausible associations of 
each class for unidentified 2FGL sources found to be positionally associated with SNRs (see 
§ 5.2.7). 

The photon flux for 1-100 GeV (-F 35 ; the subscript ij indicates the energy range as 10* - 
lO- 5 MeV) and the energy flux for 100 MeV to 100 GeV in Table H] are evaluated from the fit 
to the full band (see § 13.51) . rather than sums of band fluxes as in 1FGL. We do not present 
the integrated photon flux for 100 MeV to 100 GeV (see § 13.51) . Table [7] presents the fluxes 
in individual bands as defined in § 13.51 

Figure [18] illustrates where the different classes of sources are located in the sky. Fig- 
ure [12] shows where the bro ad classes of sources appear in the curvature - variability plane. 



This is similar to Figure 8 of lAbdo et al.l (l2010gl ) although the two indicators were improved 



Most "other" curved non-variable sources are tentatively associated to SNRs. The two "pul- 
sars" above the variability threshold are the Crab and PSR J1142+01. The Crab mixes th e 



pulsar and the nebula, and we know the variability is due to the nebula (jAbdo et al.ll201ld ). 



PSR J1142+01 is a newly discovered millisecond pulsar with no known LAT pulsations. 



4.2. Comparison with 1FGL 



The 1FGL catalog (jAbdo et al.l l2010gn lists 1451 sources detected during the first 11 
months of operation by the LAT. Associations between 2FGL and 1FGL sources are based 
on the following relation: 

A< d x = J 6* +9* (5) 



where (A) is the angular distance between the sources and d x is defined in terms of the 
semi-major axis of the x% confidence error ellipse for the position of each source, e.g., the 
95% confidence error for the automatic source association procedure (§ 15. 2p . In total, 1099 
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Note. — R.A. and Decl. are celestial coordinates in J2000 epoch, I and b are Galactic coordinates, in degrees; 6\ and O2 are the semimajor and semiminor axes of the 95%' confidence source 
location region; <j> is the position angle in degrees east of north; F35 and AF35 are photon flux 1 GeV — 100 GeV in units of 10 — 9 cm" 2 s — 1 ; 525 an d AS25 are the energy flux 100 MeV 
- 100 GeV in units of 10 — 12 erg cm -2 s — 1 ; T25 and AT25 are the photon power-law index and uncertainty for a power-law fit; Mod is the spectral model used (PL for power-law, EC for 
exponential cutoff, and LP for log parabolic); Var is the variability flag (see the text); Flags are the analysis flags (see the text); 7-ray Assoc. lists associations with other catalogs of GeV 
7-ray sources; TeV indicates an association with a point-like or small angular size TeV source (P) or extended TeV source; Class designates the astrophysical class of the associated source 
(see the text); ID or Assoc. lists the primary name of the associated source or identified counterpart; Ref. cross references LAT collaboration publications. This table is published in its 
entirety in the electronic edition of the Astrophysical Journal Supplements. A portion is shown here for guidance regarding its form and content. 



Table 5. LAT Second Catalog Description 



Column Description 



Name 2FGL JHHMM.m+DDMM[c/e] , constructed according to IAU Specifications for Nomenclature; m is decimal 

minutes of R.A.; in the name R.A. and Decl. are truncated at 0.1 decimal minutes and 1', respectively; 
c indicates that based on the region of the sky the source is considered to be potentially confused 
with Galactic diffuse emission; e indicates a source that was modeled as spatially extended (see § 13,41 1 

R.A. Right Ascension, J2000, deg, 3 decimal places 

Decl. Declination, J2000, deg, 3 decimal places 

1 Galactic Longitude, deg, 3 decimal places 
b Galactic Latitude, deg, 3 decimal places 

Q\ Semimajor radius of 95% confidence region, deg, 3 decimal places 

02 Semiminor radius of 95% confidence region, deg, 3 decimal places 

<j> Position angle of 95% confidence region, deg. East of North, decimal places 

<7 Significance derived from likelihood Test Statistic for 100 MeV-100 GeV analysis, 1 decimal place 

F35 Photon flux for 1 GeV-100 GeV, 10" 9 ph cm s , summed over 3 bands, 1 decimal place 

AF35 1 a uncertainty on F35 , same units and precision 

525 Energy flux for 100 MeV-100 GeV, 10 — 12 erg cm -2 s — 1 , from power-law fit, 1 decimal place 

AS25 1 a uncertainty on 525, same units and precision 

r Photon number power-law index, 100 MeV-100 GeV, 2 decimal places 

AT 1 a uncertainty of photon number power-law index, 100 MeV-100 GeV, 2 decimal places 

Mod. PL indicates power-law fit to the energy spectrum; LP indicates log-parabola fit to the energy spectrum; 

EC indicates power-law with exponential cutoff fit to the energy spectrum 
Var. T indicates < 1% chance of being a steady source; see note in text 

Flags See Table \3\ for definitions of the flag numbers 

7-ray Assoc. Positional associations with 0FGL, 1FGL, 3EG, EGR, or 1AGL sources 

TeV Positional association with a TeVCat source, P for angular size <20', E for extended 

Class Like 'ID' in 3EG catalog, but with more detail (see Table[6]|. Capital letters indicate firm identifications; 

lower-case letters indicate associations. 
ID or Assoc. Designator of identified or associated source 
Ref. Reference to associated paper (s) 
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Table 6. LAT 2FGL Source Classes 



Description 


Identified 


Associated 






Designator 


Number 


Designator Number 


Pulsar, identified by pulsations 


PSR 


83 






Pulsar, no pulsations seen in LAT yet 






psr 


25 


Pulsar wind nebula 


PWN 


3 


pwn 





Supernova remnant 


SNR 


6 


snr 


4 


Supernova remnant / Pulsar wind nebula 






t 


58 


Globular cluster 


GLC 





glc 


11 


High-mass binary 


HMB 


4 


hmb 





Nova 


NOV 


1 


nov 





BL Lac type of blazar 


BZB 


7 


bzb 


429 


FSRQ type of blazar 


BZQ 


17 


bzq 


353 


Non-blazar active galaxy 


AGN 


1 


agn 


10 


Radio galaxy 


RDG 


2 


rdg 


10 


Seyfert galaxy 


SEY 


1 


sey 


5 


Active galaxy of uncertain type 


AGU 





agu 


257 


Normal galaxy (or part) 


GAL 


2 


gal 


4 


Starburst galaxy 


SBG 





sbg 


4 


Class uncertain 








1 


Unassociated 








575 


Total 




127 




1746 



Note. — The designation 'f indicates potential association with SNR or PWN (see 
Table [TT}. Designations shown in capital letters are firm identifications; lower case letters 
indicate associations. In the case of AGN, many of the associations have high confidence 
l|Ackermann et al1l2011bl) . Among the pulsars, those with names beginning with LAT were 
discovered with the LAT. In the FITS version of the 2FGL catalog, the f designator is 
replaced with 'spp'; see Appendix [Cl 
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Table 7. Second LAT Catalog: Fluxes in Bands 



Name 2FGL 


100 MeV - 


300 MeV 


300 MeV - 


1 GeV 


1 GeV - 3 GeV 


3 GeV - 10 GeV 


10 GeV - 100 GeV 


Fi a 


AFi a 




F 2 a AF 2 a 




F 3 b 


AF 3 b 






AF 4 C • 


frsi 




AF 5 C 




J0000.9-0748 


1.4 


0.0 


1.5 


0.3 


0.1 


5.3 


0.5 


0.0 


2.5 


2.0 


0.0 


2.5 


0.7 


0.0 


2.0 


J0001. 7-4159 


1.5 


0.0 


2.2 


0.2 


0.0 


2.8 


0.5 


0.0 


2.6 


1.6 


0.6 


5.2 


0.9 


0.0 


0.0 


J0002. 7+6220 


1.9 


0.7 


3.5 


1.3 


0.2 


9.2 


2.6 


0.4 


8.7 


4.1 


1.1 


5.2 


0.5 


0.0 


0.0 


J0004.2+2208 


1.6 


0.0 


2.0 


0.3 


0.1 


4.8 


0.4 


0.0 


1.6 


1.3 


0.6 


3.7 


0.5 


0.0 


0.0 


J0004.7-4736 


2.2 


0.4 


5.8 


0.4 


0.1 


7.7 


0.9 


0.2 


7.9 


1.3 


0.6 


4.1 


0.6 


0.0 


0.0 


J0006.1+3821 


2.7 


0.5 


5.9 


0.5 


0.1 


7.0 


0.9 


0.2 


6.2 


2.3 


0.0 


2.8 


0.6 


0.0 


0.0 


J0007.0+7303 


17.9 


0.6 


37.6 


11.6 


0.2 


96.0 


49.9 


0.9 


122.7 


149.9 


4.2 


91.3 


12.5 


1.2 


27.1 


J0007.7+6825c 


2.8 


0.0 


0.6 


1.1 


0.2 


6.1 


0.9 


0.3 


3.2 


1.3 


0.0 


0.0 


1.0 


0.0 


1.9 


J0007.8+4713 


2.7 


0.4 


6.6 


0.6 


0.1 


8.4 


1.5 


0.2 


9.5 


4.9 


1.0 


9.3 


1.2 


0.5 


5.2 


J0008.7-2344 


0.5 


0.0 


0.0 


0.1 


0.0 


0.1 


0.6 


0.0 


3.2 


1.4 


0.0 


1.6 


1.8 


0.0 


3.8 


J0009.0+0632 


2.1 


0.0 


2.8 


0.3 


0.0 


2.4 


0.3 


0.1 


3.2 


1.4 


0.6 


3.7 


0.7 


0.0 


1.8 


J0009.1+5030 


0.9 


0.0 


0.3 


0.4 


0.1 


6.0 


1.5 


0.2 


9.6 


5.9 


1.1 


10.6 


1.9 


0.6 


7.0 


J0009.9-3206 


0.7 


0.0 


0.1 


0.3 


0.0 


3.1 


0.5 


0.1 


5.0 


1.7 


0.0 


1.6 


0.9 


0.0 


1.7 


J0010.5+6556c 


2.7 


0.0 


2.8 


1.2 


0.2 


7.0 


1.8 


0.0 


3.0 


3.2 


0.0 


1.5 


0.7 


0.0 


0.5 


J001 1.3+0054 


1.3 


0.0 


0.8 


0.3 


0.1 


4.9 


0.5 


0.1 


4.4 


1.9 


0.0 


2.0 


0.6 


0.0 


0.0 



Note. — This table is published in its entirety in the electronic edition of the Astrophysical Journal Supplements. A portion 
is shown here for guidance regarding its form and content. 

a In units of 10~ 8 photons cm -2 s~ 1 

b In units of 10 — 9 photons cm -2 s~ 1 

c In units of 10 -10 photons cm -2 s~ 1 
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Fig. 18. — Full sky map (top) and blow-up of the inner Galactic region (bottom) showing 
sources by source class (see Table [6]). Identified sources are shown with a red symbol, 
associated sources in blue. 
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Fig. 19. — Variability index (TS var in § 13. 61) plotted as a function of the curvature significance 
(Signif _Curve in § 13. 5j) for different broad classes of sources. "AGN" here means any class 
starting with "ag" or "bz" in Table |6j The horizontal dashed line is set to 41.6, above which 
sources are likely variable. The vertical dashed line is set to 4.0, above which curved spectra 
are used. 
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2FGL sources were automatically associated with entries in the 1FGL catalog. At the level 
of overlapping 95% source location confidence contours the 2FGL catalog contains 774 (out 
of 1873) new 7-ray sources, while 352 sources previously listed in 1FGL do not have a 
counterpart in the 2FGL catalog. 

The Galactic latitude distributions of the 2FGL sources, the 1FGL sources and of the 
sources in common between the 1FGL and 2FGL catalogs, shown in Figure |20| indicate both 
that most of the new 2FGL sources and most of the missing 1FGL sources are concentrated 
along the Galactic plane where the Galactic diffuse emission is most intense and improve- 
ments in the model for the diffuse emission since the 1FGL analysis would be expected to 
have the most influence (§ 12. 2j) . 

As described in § 13.31 m the 2FGL analysis the spectral fits are made using power- 
law, power-law with an exponential cutoff, or log-parabola models. Of the 1099 1FGL 
sources associated with 2FGL sources, 274 of the brightest were fitted with a curved spectral 
functional form. For each 2FGL source we also evaluated the spectral index (r) of the best 
power-law fit (§ 14. ip and this enables a comparison of the spectral characteristics of the 
1FGL and 2FGL sources. Figure [2T] shows the distributions of the spectral indices of all of 
the sources in the 1FGL and 2FGL catalogs. The two distributions are very similar, with 
an average Tifgl = 2.23 ± 0.01 and an average T 2 fgl = 2.21 ± 0.01. However, the peaks of 
the two distributions are not exactly coincident; also, the skewness of the 2FGL distribution 
is positive while it is negative for 1FGL. Figure [22] shows the distribution of the difference 
^2Fgl — Tifgl f° r the 1099 sources in common between the catalogs. The average of the 
distribution is —0.07 ± 0.01, with the 2FGL sources slightly harder than the 1FGL ones. 
This small difference in the spectral index distribution could be related to slightly difference 
uncertainties in the effective area between P7_V6 and P6_V3. 

The distributions of the source significances reported in Figure [22] show that for the 
2FGL catalog the significance peaks between 4 o and 5 a while for 1FGL the distribution 
shows a plateau between 4 a and 6 cr; this indicates that 2FGL is more complete than 1FGL. 
Also, the distribution of the significance of the sources that are in common between 1FGL 
and 2FGL shows that most of the 1FGL sources that were not recovered in the 2FGL catalog 
had significance less than 7 a. In the remainder of this section we describe the variety of 
reasons that the additional 352 1FGL sources do not appear in the 2FGL catalog. 

Table M lists 347 of the 1FGL sources that do not have a corresponding source in 
2FGL. The five other 1FGL sources that do not appear in 2FGL were not included in 
the table because they were already replaced by an extended source template in the 2FGL 
analysis. These sources are: 1FGL J0523.3-6855 (2FGL J0526.6-6825e, LMC); 1FGL 
J1801.3-2322c (0FGL J1801. 6-2327, 2FGL J1801.3-2326e, W28); 1FGL J1805.2-2137c 



-59 - 



180 
160 
140 
120 
100 
80 
60 
40 
20 




1FGL-2FGL common sources 

2FGL 

1FGL 




-1.0 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1.0 

sin(GLAT) 



Fig. 20. — Distributions of the Galactic latitude of the 1FGL and 2FGL sources and of the 
sources in common between the 1FGL and 2FGL catalogs. 



(2FGL J1805.6-2136e, W30); 1FGL J1856. 1+0122 (2FGL J1855.9+0121e, G034.7-00.4, 
W44); 1FGL J1922.9+1411 (2FGL J1923.2+1408e, G049.2-00.7, W51C). 

Some 1FGL sources near extended 2FGL sources remain in Table El An additional 
four 1FGL sources, 1FGL J0459.7-6921, 1FGL J0518.6-7222, 1FGL J0531. 3-6716 and 
1FGL J0538.9— 6914, were also found in the LMC field, now replaced by an extended source 
in the 2FGL catalog analysis (J0526.6— 6825e). Furthermore, the 4 1FGL sources 1FGL 
J2046.4+3041, 1FGL J2049. 1+3142, 1FGL J2055.2+3144, 1FGL J2057.4+3057 distributed 
along the Cygnus Loop (G74. 0—8.5), one of the most famous and well-studied SNRs, were 
replaced by an extended source template in the 2FGL analysis (2FGL J2051.0+3040e), and 
so are not confirmed in the 2FGL catalog. The extended source 2FGL J1824.5— 1351e (HESS 
J1825-137) replaces two 1FGL sources: 1FGL J1821.1-1425c and 1FGL J1825. 7-1410. 

About 250 of the 347 sources are located on the Galactic plane or in other regions of 
bright, structured diffuse emission (see Fig. |2"4]) . Of these, 88 have the 'c' designation in the 
1FGL name, which indicates that these sources were already recognized as possible spurious 
detections. Another 21 1FGL so urces were flagged according to the definitions reported 
in Table 4 of lAbdo et al.l (I2010gl ). These sources were also already noted as problematic. 
In the 1FGL catalog only 67 of the 347 sources have an association with a possible coun- 
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Fig. 21. — Distributions of the spectral index for the 1FGL (1451 sources, dashed line) and 
for the 2FGL (1873 sources, solid line) catalogs. 



terpart, mostly AG N, while a nother 10 sources w ere associated with already known OFGL 
jAbdo et al.ll2009dh or 1AGL Jpittori et al-lbopgh 7-ray sources. 



In addition to the introduction of spatially extended sources in 2FGL, there are many 
possible causes for 1FGL sources to be absent from the 2FGL list. Among these are variabil- 
ity; different event selection used for the analysis (Pass 6 for 1FGL and Pass 7 for 2FGL); 
different IRFs; different Galactic diffuse emission models; different analysis procedures (un- 
binned likelihood analysis for 1FGL and binned likelihood analysis for 2FGL); statistical 
threshold effects; and 1FGL sources resolved into two or more 2FGL sources. In the last 
columns of Table [8] we assigned to each source one or more flags corresponding to a possible 
cause. In many cases, no one reason can be singled out. 

The numbers of associated sources between the 1FGL and 2FGL catalogs does depend 
on the criterion used to define spatial coincidence (Eon. 5). The number of 2FGL - 1FGL 
associated sources increases to 1151 if we use A < g? 99 .cl|. The 52 additional associations (see 
Table Eland see Figure 123]) . represent about the 5% of the 1451 1FGL sources, as expected 



Assuming a Rayleigh distribution for the source angular separations, c/99.9 is evaluated using ^99. 9 = 
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Fig. 22. — Distribution of the difference T 2 fgl — ^ifgl for the 1099 sources in common 
between the 1FGL and 2FGL catalogs. 



when passing from dg§ to dgg.g. Also, in 2FGL we used a better in-flight representation 
of the PSF that is broader than the PSF used in 1FGL at energies E> lGeV where, in 
general, most of the sources are detected. Furthermore, the new and improved model of 
the Galactic diffuse emission used to build the 2FGL catalog together with the expected 
increase of the signal-to-noise ratio due to the use of 24 months data, allowed us to obtain 
a better localization of the sources at positions that might be outside the 95% confidence 
error regions previously reported in 1FGL. Indeed, most of the 52 additional associations 
concern sources located along the Galactic plane and in regions like Orion and Ophiuchus, 
while only about 10 were associated in regions with low diffuse emission. 

Also, in the 1FGL catalog the positions of sources associated with the LAT-detected 
pulsars and X-ray binaries are the high-precision positions of the identified sources. (These 
sources can be easily recognized because they have null values in the localization parameters 
reported in the 1FGL catalog). Not all of these associations appear in the 2FGL catalog 
because they cannot be associated using dg 5 , but some are listed in Table [8] because they 
can be associated using dgg,g. These sources are: 1FGL J2032. 4+4057 (Cyg X— 3); 1FGL 
J1836.2+5925 (LAT PSR J1836+5925); 1FGL J1124.6-5916 (PSR J1124-5916). However, 



1.52 6, 
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Fig. 23. — Distributions of the significances of the sources in the 1FGL and 2FGL catalogs 
and of the sources in common between the catalogs. The distribution of the significance for 
the 1FGL-2FGL common sources is based on the values reported in the 1FGL catalog. 



1FGL J1741.8-2101 (LAT PSR J1741-2054), 1FGL J1614.0-2230 (PSR J1614-2230) and 
1FGL J1747.2-2958 (PSR J1747-2958) are still not associated, and for these sources we 
report the nearest 2FGL source (see, e.g., Fig. [26]). The last missing source in this category 
is 1FGL J1023.0-5746 (LAT PSR J1023-5746). It was resolved into two 2FGL sources, 
2FGL J1022.7-5741 and 2FGL J1023.5-5749. Although both are located very close to the 
pulsar position, they cannot be formally associated using d 99 , 9 . 

Several other 1FGL sources were also split into more than one candidate source seed ('S', 
in the Flags column of Table [8]). In some cases only one of the two seeds reached a TS>25 
and so was included in the 2FGL list (see Figure 1271) . Another example of splitting i s 1FGL 



J1642. 5+3947, that was tentatively associated with the blazar 3C 345 in lAbdo et al.l (1LAC 



2010a ) paper. This source has no 2FGL counterpart, because it is now resolved into two 
sources: 2FGL J1642.9+3949 associated with 3C 345 and 2FGL J1640.7+3945 associated 
with NRAO 512. Other 1FGL sources have overlapping #99.9 source location uncertainty 
regions with one or more 2FGL sources or seeds and have the 'O' flag in Table [H] (see Fig. 



Another major reason for sources to disappear between 1FGL and 2FGL is a change in 
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Fig. 24. — All-sky map for energies >lGeV indicating the positions of the 1FGL sources 
that are not in the 2FGL catalog (green circles). The red crosses indicate the sources having 
the flag £ NC in Table E) 

the calculated significance. As described in the § I3.1.3} the 2FGL catalog was built starting 
from 3499 seeds with TS > 10 in the pointlike analysis. The final gtlike analysis, which 
did not change the positions of the seeds, resulted in the 1873 sources with TS > 25 that 
make up the 2FGL catalog. Among the other seeds that did not reach the threshold, 104 
can be associated with 1FGL sources (using A < ^99.9). These sources, marked with a 'C in 
the flags column of the table, can be considered to be confirmed sources whose significance 
dropped below the threshold, either as a result of time variability, change in the diffuse 
model, or the shift from unbinned to binned likelihood in the catalog analysis procedure. 

In order to quantify the effect of changing gtlike from unbinned to binned mode, we 
performed a new binned analysis of the original 11-month data set, using the P6V3 Diffuse 
IRFs and the same Galactic diffuse emission model as used for the 1FGL analysis. The 
analys is also started usin g the same 1499 seeds that were used as input to the 1FGL run 



(see, lAbdo et al.l l2010gl ). This analysis confirmed with TS>25 1138 sources of the 1451 
sources that were in the 1FGL catalog. Among these confirmed sources are 168 1FGL 
sources that are not present in 2FGL, but were still detected at TS> 25 using the binned 
analysis for the 11-month data set. In Table |8]the sources confirmed by the binned analysis 
but not included among the 3499 seeds have the flag 'BC\ In the 1FGL catalog, only 5 of the 
'BC sources were found to be variable with probabilities p > 90% (see the 'Var' column in 




Fig. 25. — A typical example of a 1FGL source associated with a 2FGL source using d$g,$. 
The E > 1 GeV counts map (1 pixel = 0?2) was smoothed using a gaussian kernel (a = 3 
pixels). The ellipses represent the 99.9% confidence error regions. 



Table [8]). Since the shift from unbinned to binned analysis has been excluded for 
these, their disappearance must be attributed to time variability or, more likely, to change 
in the diffuse emission model. 

The 102 sources that were no longer detected in either the binned likelihood re-analysis of 
the 11-month dataset nor in any of the other all-sky analyses performed using data collected 
between 11 and 24 months are considered not confirmed 1FGL sources ('NC in the Flags 
column of Table ED- Among these sources 9 were bright during just the first months of the 
mission and are reported with the flag 'V (variable) in Table |HJ They are all associated with 
AGN, mostly blazars. An exa mple is 1FGL J1122 9— 6415, associated with PMN J1123-6417 



and included in the OFGL list ( lAbdo et al.ll2009dl ). that, after a flare in 2008 September was 
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Fig. 26. — A typical example of a 1FGL and a 2FGL source having an angular separation 
greater than dgg.g but less than 1°. In this particular case there is also a seed (candidate 
source considered in the 2FGL analysis) separated by less than 1° from the 1FGL source. 
The E > 1 GeV counts map (1 pixel = 0?2) was smoothed using a gaussian kernel (a = 3 
pixels). The ellipses represent the 99.9% confidence error regions. 



not significantly detected again until 2011 May ( ID ' Ammandol l201lh . 



Four 'NC 1FGL sources have angular distances less than 1° from the ecliptic ('Sun' in 
the Flags column of Table E]). Their light curves, which are similar to that shown in Figure 
[TBI show significant detections of the sources only during the passage of the Sun. For the 
other 89 'NC sources their non-associations with the 2FGL sources can be ascribed to a 
combination of different effects that cannot be easily disentangled. 



Most of these sources are located close to regions of enhanced diffuse emission (see 
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Fig. 27. — A typical example of a 1FGL source that was split in two 2FGL sources. The E 
> 1 GeV counts map (1 pixel = 0?2) was smoothed using a gaussian kernel (a = 3 pixels). 
The ellipses represent the 99.9% confidence error regions. 

Figure EIJ) and about 20 of them were already flagged as sources influenced by the diffuse 
emission in the 1FGL catalog. Also, the fact that these sources were not confirmed in the 
binned analysis of the 11-month data can be related to statistical fluctuation in the number 
of the sources detected close to the significance threshold. Figure [29] shows the distribution 
of the source significances, as reported in 1FGL, for the 89 non-confirmed 1FGL sources and 
for the 1099 1FGL sources present in the 2FGL catalog. Most of the non-confirmed sources 
have significances less than 6 a, which is very close to the threshold ( ~ 4 a) adopted in 1FGL 
and 2FGL catalogs. These sources are intrinsically faint and for several of them the energy 
flux reported in the 1FGL catalog is just an upper limit. Furthermore the significance values 
returned by the unbinned likelihood analysis by definition should be intrinsically higher than 
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Fig. 28. — A typical example of a 1FGL source and a 2FGL source having overlapped 
99.9% confidence error regions. In this particular case there is also a seed (candidate source 
considered in the 2FGL analysis) very close to the 1FGL source and we cannot exclude the 
possibility that the 1FGL source was split into two seeds. The E > 1 GeV counts map (I 
pixel = 0?2) was smoothed using a gaussian kernel (er = 3 pixels). The ellipses represent the 
99.9% confidence error regions. 

those returned by the binned analysis. Thus, most of the 89 sources were above the threshold 
in the original unbinned 1FGL analysis, but not in the binned analysis. 



Table 8. List of 1FGL sources not in the 2FGL catalog 
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a All the values reported in these columns are from the 1FGL catalog l lAbdo et al.ll2010gh . 
b Name of the 2FGL source associated with the 1FGL one using dgg.g. 

c Closest 2FGL source having a distance dgg.g < A < 1° from the position of the 1FGL source. The 2FGL name is also reported if the 1FGL source and one or more 
seeds have overlapping S99 9 error regions but cannot be associated with any seed on the basis of the criterium A < cigg.g. 

d The angular separation (A) between the 1FGL source and the 2FGL sources associated using dgg.g or the closest 2FGL source. 

°T = The 1FGL source and one of the 2FGL list of initial seeds have an angular separation < dgg.g. 

f C= Confirmed 1FGL sources. LMC, Orion, Carina and Ophiuchus indicate that the source is in a region of the sky with high diffuse emission and high density of close 
sources; NC = not confirmed 1FGL sources (see text); BC = 1FGL sources confirmed by the 11-m binned likelihood analysis; S = the 1FGL source was split/resolved in 
one or more seeds; O = overlapping Sgg.g error regions with one or more seeds; V = variable source visible only in the first 11 months; Sun = the source was detected when 
the Sun was at an angular distance < 1° and the light curve show just a flare in the time bin relative to the passage of the Sun close to the position of the source. This 
table is published in its entirety in the electronic edition of the Astrophysical Journal Supplements. A portion is shown here for guidance regarding its form and content. 
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Fig. 29. — Distribution of the significances of the unconfirmed 1FGL sources and of the 
1FGL sources associated with sources in the 2FGL catalog. 

5. Source Association and Identification 



5.1. Firm Identifications 



As with the LAT Bright Source List f lAbdo et alj|2009dj ) and 1FGL catalog flAbdo et al. 



2010gj), we retain the distinction between associations and firm identifications. Although 
many associations, particularly those for AGN, have very high probability of being true, a 
firm identification, shown in the catalog by capitals in the Class column in Table El is based 
on one of three criteria: 



1. Periodic Variability. Pulsars are the larger class in this category. All PSR labels indi- 
cate that pulsed 7 rays have been seen from the source with a probability of the period- 
icity occurring by chance of less than 10~ 6 . A similar chance probability requirement 
applies to the other set of periodic sources, t he high-mass binari es (HMB). Four of these 



are in cluded in the catalog: LS I +61 303 flAbdo et al.ll2009d ). L S 5039 flAbdo et al. 



2009eJ), Cygnus X-3 flAbdo et al.ll2009ff ) . and 1FGL J1018.6-5856 f lCorbet et al.ll201lh . 
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2. Spatial Morphology. Spatially extended sources whose morphology can be related to 
extent seen at other wavelengths include SNR, PWNe, and galaxies, as described in 
§ 13.41 The Centaurus A lobes and core are both marked as identified, because they are 
part of the same extended source, although the core itself does not show spatial extent. 
As noted in § 3.8, additional extended sources are being found but are not listed in 
the catalog as firm identifications, because they were analyzed as point sources for this 
work. Although individual molecular clouds could in principle be included in this list, 
the catalog construction incorporates most known clouds into the diffuse model, and 
so no sources of this type are identified in the catalog. 

3. Correlated Variability. Variable sources, primarily AGN, whose 7-ray variations can be 
matched to variability seen at one or more other wavelengths, are considered to be firm 
identifications. Although some cases are well documented, such correlated variability 
is not always easily defined. We conservatively require data in more than two energy 
bands for comparison. Finding a blazar to have a high X-ray flux at the same time as a 
7-ray flare, for example, does not qualify if there is no long-term history for the X-ray 
emission. We include those sources whose variability properties are documented either 
in papers or with Astronomer's Telegrams. This list does not represent the result of a 
systematic study. Ongoing work will undoubtedl y enlarge this list. The one Galactic 



source identified in this way is nova V407 Cygni ( lAbdo et al.ll2010rj ). 



We include one exception to these rules. The Crab PWN is listed as a firm identification 
even though it does not meet any of these criteria. The well-defined energy spectrum, distinct 
from the Crab pulsar spectrum and mat ching spectr a seen at both lower and higher energies 



provides a unique form of identification lAbdo et al.l ( 12010jl ) . 



In total, we firmly identify 127 out of the 1873 2FGL sources. Among those, 83 are 
pulsars, 28 are AGN, 6 are SNR, 4 are HMB, 3 are PWN, 2 are normal galaxies, and one is 
a nova (Table [6]). 



5.2. Automated Source Associations 



Our approach for automated source associat ion closely follows t hat used for the 1FGL 
catalog, and details of the method are provided in lAbdo et al.l ( l2010gT ) . In summary, we use a 
Bayesian approach that trades the positional coincidence of possible counterparts with 2FGL 
sources against the expected number of chance coincidences to estimate the probability that 
a specific counterpart association is indeed real (i.e., a physical association). As for 1FGL, 
we retain counterparts as associations if they reach a posterior probability of at least 80%. 
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We apply this method to a set of counterpart catalogs for which we calibrate the prior 
source association probabilities using Monte Carlo simulations of fake 2FGL catalogs. In 
comparison to 1FGL, for which we made 100 independent simulations for each catalog, we 
adapted the number of simulations (between 100 and 1000) so that the relative accuracy 
in the expected false association rate is determined to better than 5% for each catalog. 
This improved the precision of our probability computations for catalogs that have only few 
associations with 2FGL sources. The prior probabilities adopted for each catalog are listed 
in Table El 

Another improvement with respect to 1FGL concerns the estimation of the local coun- 
terpart densities pt- For 1FGL we estimated these densities from the number of objects in 
the counterpart catalog within a radius of 4° around the location of the 1FGL source of inter- 
est. For counterpart catalogs containing strong density variations on smaller scales (e.g., O 
stars, WRs and LBV stars) this choice led to an underestimate of the actual source densities 
in these regions, which in turn resul ted in overestimations of the association probabilities 



(see discussion in lAbdo et al.1 l2010gl ) . For 2FGL we estimate the source densities in each 
counterpart catalog using an all-sky map which we implemented as a HEALPix grid with 
resolution N side = 512, corresponding to an angular resolution of about 6', with the objects 
of each counterpart catalog binned in this grid. We removed sparseness of the binning and 
attenuated the statistical fluctuations by applying a spherical Gaussian smoothing kernel 
with width adjusted adaptively so that at least 3 sources contributed to the density estimate 
at each grid location. 

For certain counterpart catalogs the Bayesian method could not be applied since either 
(1) the location uncertainty of the counterpart is larger than the location uncertainty of the 
2FGL source (these catalogs are indicated by * in Table [9]), or (2) the counterpart is an 
extended source (these catalogs are indicated by t in Table [9]). In the first case, we consider 
as potential associations all objects for which the separation from the 2FGL source is less 
than the quadratic sum of the 95% confidence error radii. (For elliptical error regions we take 
the semimajor axis as the error radius.) In the second case, we assume that the counterparts 
have circular extensions and consider all objects as associations for which the extension circle 
overlaps with the 95% confidence error radius of the 2FGL source, with the semimajor axis 
of the 2FGL source location ellipse again taken as the error radius. 

The list of catalogs used in the automatic association is summarized in Tabled organized 
into four categories: (1) catalogs of known or plausible 7-ray-emitting source classes, (2) 
catalogs of surveys at other frequencies, (3) catalogs of GeV sources, and (4) catalogs of 
identified 7-ray sources. The first category allows us to assign 2FGL sources to object 
classes, while the second category reveals multiwavelength counterparts that may suggest 
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the possible nature of the associated 2FGL source. The third category allows assessment 
of former GeV detections of 2FGL sources, and the fourth category keeps track of all firm 
identifications (cf. 



15. ip . For this last category we claim associations based on the spatial 
overlap of the true counterpart position with the 2FGL 99.9% confidence error ellipse. 

With respect to 1FGL, we updated all catalogs for which more com prehensive compila- 



tions became available. We no w use the 13th editio n of the Veron catalog (IVeron-Cettv fc Veron 
201oh . version 20 of BZCAT@ jMassaro et alil2009h . version 1.40 of ATNF jManchester et al 



20051 ) that we augmented with 158 recently detected pulsars that are not yet in the A TNF 
database, the 2010 December revision of the Globular Cluster database (jHarria 119961 ) that 
we augmented w ith 3 recently detected clusters, version 3.1 of the Open Cluster catalog 
( Dias et al.l 12002 ). the 2010 December 5 version of the VLBA Calibrator Source Lisfl and 
the most recent version of the TeVCat c atalog. We also ad ded new counterpart catalogs: 
the Australia Te lescope 20 GHz Sur vey (IMurphy et al.l 120101 ) and the IRAS Revised Bright 
Galaxy Sample (jSanders et al.l 120031 ) . from which we selected all sources with 100 /im fluxes 
brighter than 50 Jy. The latter catalog replaces the starburst catalog used for 1FGL. 

Following the philosophy for 1FGL, we split our pulsar catalog into normal pulsars and 
millisecond pulsars (MSPs) by requiring log P + 19.5 + 2.5 x logP < for the latter. Because 
globular clusters are classified by a separate catalog and the LAT is unable to spatially resolve 
individual MSPs in globular clusters, we removed all globular cluster MSPs from the pulsar 
catalog. We furthermore collect normal pulsars with E/d 2 > 5 x 10 32 erg kpc -2 s -1 into 
a separate counterpart catalog to specifically select energetic and nearby pulsars that are 
more likely potential 7-ray sources. The value separating these classes corresponds to the 
lowest E/d 2 found among all LAT identi fied pulsars. We also split-off point-like supernova 
remnants (SNRs) from the Green catalog (lGreenll2009i ) by selecting all objects with diameters 
< 20'. In parallel, we use the full Green catalog for finding matches with potentially extended 
SNRs. Furthermore we divided the TeVCat catalog into point-like and extended sources by 
selecting for the latter all sources with extension radius > 0. 

We also searche d for associations using the Atlas of Radio/X-ray associations to optical 
objects (lFleschll2010l ) from which we selected those objects th at have stellar, ra dio, and X-ray 



associati ons (C1=SRX), t he Planck Early Release Catalogs (lAde et al. 



catalog (Bird et al.ll2010h . and the Swift-BAT 58- Month Survey (IBaumgartner et all 120101 ) 



201 lah. the 4th IBI S 



9 http://www. asdc.asi.it/bzcat/ 

10 The VLBA Calibrator Source 
http : / / www . vlba. nrao . edu/astro/calib / vlbaCalib . txt . 

11 http:/ /tevcat. uchicago.edu/ 
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yet as these did not reveal any new reliable and plausible counterpart that has not already 
been found in one of the other catalogs, we did not include these catalogs in our final results. 



5.2.1. Automated association summary 

The results of the automated association procedure for each of the external catalogs are 
summarized in Table [9) For each catalog we quote the name (Column 1), the number of 
objects in the catalog (Column 2), the prior probability assigned by our calibration proce- 
dure (Column 3), and the number iV ass of associations that have been found between 2FGL 
sources and counterpart objects (Column 4). Note that a given 2FGL source may have coun- 
terparts in multiple catalogs, and a given object in a counterpart catalog may have multiple 
associated 2FGL sources (which may arise if the object is spatially extended or if it has a 
large location uncertainty). Consequently, the sum of the iV ass column considerably exceeds 
the total number of associated 2FGL sources. Using the posterior probabilities that we 
derive by the Bayesian method for all associations % in a counterpart catalog k, we compute 
the expected number of false associations using N {alse = Yup- 0- ~~ ^»*) (Column 5). To vali- 
date that these estimates are accurate (and thus that our prior probability calibration was 
precise) we alternatively estimate the number of false associati ons {.\'r :< u,.) u s ing Mo nte Carlo 



simulations of 100 fake 2FGL catalogs (Column 6); we refer to lAbdo et al.l (l2010gl ) for a de- 
tailed description of the simulation procedure. For all catalogs we find iVf a i se ~ (Nf&\ S e) which 
confirms that the posterior probabilities computed by the automatic association procedure 
are accurate. 

In total we find that 1141 of the 1873 sources in the 2FGL catalog (61%) have been 
associated with a least one non-GeV 7-ray counterpart by the automated procedure. Among 
those, 123 sources (11%) are firmly identified objects, 790 (69%) are associated with at least 
one object of known type, and 228 (20%) have counterparts only in the mult i- wavelength 
catalogs. For the remaining 732 sources in the 2FGL catalog that have no non-GeV 7-ray 
counterpart, 322 sources (44%) are associated with former GeV detections, and 410 sources 
(56%) are new GeV sources. 

Among the 2FGL sources that are not firmly identified, 940 (92%) have been associated 
using the Bayesian method at the 80% confidence level, while 78 (8%) have been associated 
based on overlap of the error regions or source extents and have lower confidence (catalogs 
based on spatial overlap are indicated by ' in Table EJ). From simulations we expect that 43 of 
the 940 sources (5%) that were associated with the Bayesian method are chance coincidences. 
Among the 78 sources that were associated based on overlap, the expected number of chance 
coincidences amounts to 55 (71%), demonstrating that these associations are considerably 
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Table 9. Catalogs used for the automatic source association 



Name 


Objects 


P prior 


Nass 




(Afalsc) Ref. 


High E/d 2 pulsars 


213 


0.037 


29 


0.9 


1.0 


2 


Other normal pulsars 


1657 


0.011 


12 


0.6 


0.7 


2 


Millisecond pulsars 


137 


0.014 


45 


0.3 


0.4 


2 


Pulsar wind nebulae 


69 


0.009 


25 


0.5 


0.6 


1 


High-mass X-ray binaries 


114 


0.003 


2 


0.1 


0.2 


3 


Low-mass X-ray binaries 


187 


0.007 


3 


0.3 


0.3 


1 


Point-like SNR 


157 


0.019 


6 


0.7 


0.3 


5 


Extended SNPj 


274 


n.a. 


92 


n.a. 


39.7 


5 


O stars 


378 


0.005 


1 


0.2 


0.2 


6 


WR stars 


226 


0.005 








0.2 


7 


LBV stars 


35 


0.001 


1 


< 0.1 


0.2 


8 


Open clusters 


2140 


0.005 








0.2 


9 


Globular clusters 


160 


0.028 


11 


0.5 


0.6 


10 


Dwarf galaxies t 


14 


n.a. 


7 


n.a. 


3.4 


1 


Nearby galaxies 


276 


0.014 


5 


0.4 


0.4 


11 


IRAS bright galaxies 


82 


0.021 


6 


0.2 


0.2 


12 


BZCAT (Blazars) 


3060 


0.341 


691 


7.4 


6.9 


13 


BL Lac 


1371 


0.170 


278 


2.8 




14 


AGN 


10066 


0.009 


8 


0.3 


0.4 


11 






u. lyo 


1 (17 


0. / 


6.7 


14 


Seyfert galaxies 


27651 


0.028 


29 


2.0 


1.9 


14 


Radio loud Seyfert galaxies 


29 


0.001 


4 


< 0.1 


< 0.1 


1 


CGRaBS 


1625 


0.258 


352 


3.8 


4.1 


15 


CRATES 


11499 


0.341 


634 


17.7 


17.8 


16 


VLBA Calibrator Source List 


5776 


0.258 


623 


11.8 


12.0 


17 


ATCA 20 GHz southern sky survey 


5890 


0.296 


335 


10.3 


10.6 


18 


TeV point-like source catalog* 


61 


n.a. 


47 


n.a. 


0.6 


19 


TeV extended source catalogt 


57 


n.a. 


48 


n.a. 


20.1 


19 


1st AGILE catalog* 


47 


n.a. 


57 


n.a. 


21.1 


20 


3rd EGRET catalog* 


271 


n.a. 


116 


n.a. 


31.0 


21 


EGR catalog* 


189 


n.a. 


69 


n.a. 


11.4 


22 


OFGL list* 


205 


n.a. 


185 


n.a. 


5.1 


23 


1FGL catalog* 


1451 


n.a. 


1099 


n.a. 


18.1 


24 


LAT pulsars 


87 


n.a. 


80 


n.a. 


1.4 


1 


LAT identified 


44 


n.a. 


43 


n.a. 


0.7 


1 



References. - 1 Collaboration internal; ^Manchester et al.l J2005h; j Liu et al.l 
&200d); ^Liu et alJ J2007I); ^Greenl j200Sh ; iMai'z-Apellaniz et al] 



._ der Huchtl feOOlh; ^Clark et al 
il996h: ^Schmidt et all <1993h 
HVeron-Cettv fc Veronl j201Ch ; 



.20051); iDias et al 
^Sanders et all 200; ); 
toalev et ail hooj ); 



l 

17 http: //www. vlba.nrao.edu/astro /calib/vlbaCa li b.txt; 
1 9 http://tevcat.uchicago . edu/; 2c |Pittori et al l J20091 ) 



s |MassaroetalJ 
iHealev et ajj 



u/. 

^Casandiian fc Grenierl j200Sh : 2 ^Abdo et al.l i2009dh : 



"Murphy et al. 
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less reliable. Due to this large false positive rate, we do not claim any associations based on 
overlap in our final catalog. We record, however, any spatial overlap with a TeV source in 
the FITS file version of the catalog, and use a special flag in our catalog (TEVCAT_FLAG), 
distinguishing point-like (P) from extended (E) TeV counterparts (see Appendix C). We 
furthermore list all unidentified 2FGL sources that are spatially overlapping with SNRs in 
Table [TTJ Finally, 2FGL sources spatially overlapping with the LMC that are not associated 
with any object in one of the other counterpart catalogs are indicated as LMC field. 



5.2.2. Active Galactic Nuclei associations 



Active Galactic Nuclei (AGN), and in particular blazars, are the most prominent class 
of associated sources in 2FGL. In total, our automatic association procedure finds 917 2FGL 
sources that are associated with AGN, of which 894 are blazars, 9 are radio galaxies, 5 
are Seyfert galaxies, and 9 are other AGN. Among the 5 Seyfert galaxies, 4 are narrow-line 



Seyfer t 1 galaxies that have been established as a new class of 7-ray active AGN (lAbdo et al. 



2009ef ). The 5th object is NGC 6814, which is associated with 2FGL J1942.5-1024. Note, 
however, that we expect up to ~ 2 false positives among the Seyfert galaxy associations 
(cf. Table [9]), hence we cannot draw any firm conclusions about the possibility that normal 
Seyfert galaxies are indeed GeV 7-ray sources based on this single association. 

AGN observed by the LAT are also sources of radio (and X-ray) emission, and we find 
a clear trend that AGN associated with 2FGL sources have larger radio fluxes than the 
average object in the counterpart catalogs. Th is trend, which was exploited already for the 
association of blazars in the EGRET catalog (ISowards-Emmerd et al.l 120031 ). is illustrated 
in Figure [301 where we compare the distribution of the 8.4 GHz radio fluxes of all sources 
in the CRATES catalog to that for objects associated with 2FGL sources. Obviously, the 
average radio flux of CRATES sources associated with 2FGL sources is about one order of 
magnitude larger than the overall average for the CRATES catalog. Similar differences are 
observed for other radio catalogs. 

In our dedicated effort for studying the AGN population in the 2FGL catalog, which 
we publish in an accompanying paper (2LAC; lAckermann et al.ll2011bl ). we make use of this 
property to enhance the sample of associated 2FGL sources. Briefly, instead of including 
all objects from the counterpart catalog in the estimation of the local counterpart densities 
Pk, we count only those objects with radio (or X-ray) flux equal or larger than the flux S 
of the counterpart under consideration, i.e., pk(> S). Using this procedure, the chance co- 
incidence probabilities are consider ably reduced, and consequentl y, the posterior association 
probabilities are increased (see also ISowards-Emmerd et al.ll2003l ). We apply this procedure 
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to a number of fairly uniform surveys of radio sources (CRATES, NVSS, SUMSS, PMN, 
ATCA 20 GHz, FRBA, GAPS, CLASS and VCS) and to the ROSAT All-Sky Survey of 
X-ray sources (RASS), for which we assume the counterpart density Pk(> S) to be position 
independent. In this case pk{> S) is then determined from the logiV — log S distribution of 
objects in the catalog divided by the survey area, where N is the total number of sources 
with flux > S. 

The 2LAC association procedure increases the number of AGN associations by 173, 
resulting in a total of 1090 2FGL sources that we associate with known AGN; note that the 
2LAC catalog lists two associated AGNs for the 27 2FGL sources for which more than one 
plausible association was found. The total number of 2FGL sources associated with a least 
one non-GeV 7-ray counterpart is thus 1314 (70% of all 2FGL sources). Among the AGN as- 
sociations we find 1064 blazars, of which 432 are BL Lac (+38 with respect to the automatic 
association procedure), 370 are FSRQ (+24), and 262 are of unknown type (+108). The pro- 
cedure also reveals 2 additional radio galaxies (For A associated with 2FGL J0322.4— 3717 
and PKS 0943—76 associated with 2FGL J0942.8— 7558), and one additional Seyfert galaxy 
(ESO 323-77 associated with 2FGL J1306.9-4028). For the final AGN associations pre- 
sented in the 2FGL catalog, we adopt the results of the 2LAC procedure combined with the 
results of the automatic association pipeline (see also Table |6]). 



Comparing to 1FGL (lAbdo et al.ll2010gj ). where out of 1451 sources 573 (40%) were 



found to be associated with blazars, 802 (43%) out of 1873 sources are associated with blazars 
in 2FGL, a relative increase which is readily explained by the particular effort that ha s been 



undertaken to maximize the number of blazar associations (lAckermann et al.ll2011bl ). Ne- 
glecting the 2LAC blazar associations, the fraction of 2FGL sources associated with blazars 
would have been 40%, identical to what was found for 1FGL. On the other hand, the pro- 
portion of active galaxies of uncertain type (designated by 'agu' in Table 0J has increased 
considerably: while 92 (6%) 1FGL sources were classified 'agu', 262 (14%) 2FGL sources are 
now in this category, more than doubling the proportion of this source class. This increase 
can be explained by the extensive use of radio and X-ray surveys in the 2LAC association pro- 
cedure that provides a greater number of blazar candidates that deserve dedicated follow-up 
observations to assess their natures. We also note that in the 2FGL catalog we have two new 
extragalactic source classes with respect to 1FGL: radio galaxies ('rdg') and Seyfert galaxies 
('sey'). Both were counted in the 'non-blazar active galaxy' class (designated by 'agn') in 
1FGL, and 28 (1.9%) 1FGL sources were associated with that class. Adding the 'rgd' and 
'sey' designators to the 'agn' for 2FGL amounts to 27 (1.4%) associations, a number that is 
comparable to that found for 1FGL. 
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5.2.3. Normal Galaxies 



Normal galaxies are now established as a class of high-energy 7-ray emitters f lAbdo et al. 

2010i[ ). and we associate 7 2FGL sources with such objects. Of those, we consider the 



Small Magellanic Cloud (2FGL J0059.0-7242) and the Large Magellanic Cloud (LMC, 
2FGL J0526.6— 6825) as identified owing to their spatial extensions in the LAT data. From 
the remaining five, 4 are classified as starburst galaxies: M82 (2FGL 0955.9+6936), NGC 253 
(2FGL J0047.0-2516), NGC 4945 (2FGL J1305.8-4925), and NGC 1068 (2FGL J0242. 5+0006). 
The fifth is the Andromeda galaxy M31 (2FGL J0042.5+4114). 

Except for M31, all of the associated 2FGL 7-ray sources in this class were already 
present in 1FGL, yet the two starburst galaxies NGC 4945 and NGC 1068 were not associated 
as su ch as they were no t included in our very limited counterpart catalog used at that 
time (lAbdo et al.l l2010gf ). For 2FGL, we included a catalog of infrared bright galaxies in 
the automatic association procedure (see § 15. 2p because starburst galaxies are prominent 
emitters in this waveband. Furthermore, we have found that th e 7-ray fluxes of Lo cal Group 
and starburst galaxies correlate well with star formation rates ( lAbdo et al.ll2010il ). which in 
turn correlate with infrared luminosity. Hence by selecting infrared bright galaxies from the 
IRAS Revised Bright Galaxy Sample (ISanders et al.l 120031 ) we have added a catalog to our 
procedure that contains normal galaxies that are potential 7-ray emitters. 

Three 2FGL sources lie within the extended-source template for the LMC (2FGL J0451.8— 7011, 
2FGL J0455.8-6920, and 2FGL J0533.3-6651; see §E3D- Their physical association with 
the LMC is not certain, but they are classified here as being part of the LMC and as men- 
tioned in § 15.2. II are indicated in Table H] as belonging to LMC field. 



5.2.4- Pulsars 

As of this writing, 87 pulsars have been firmly identified by the LAT through the detec- 
tion of 7-ray pulsations. Four of these pulsars did not pass TS > 25 in the catalog analysis, 
and therefore they were excluded from the 2FGL catalog. These pulsars are PSR J1513— 5908 
(aka PSR B1509-58), PSR J1531-5610, PSR J1801-2451, and PSR J1939+2134. Of the 
remaining 83, 80 were formally associated by the automatic association procedure. The re- 
maining 3 are found to be close to 2FGL sources, but their angular separation A from these 
sources exceeds their effective 99.9% location error radius #99.90 We find: 



12 The effective error radius is the size of the error ellipse at the position angle toward the counterpart. 
We estimate the 99.9% confidence radius by multiplying the 95% confidence radius by 1.52. 
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PSR J1023-5746 near 2FGL J1022. 7-5741 (6 99 . 9 = 4.2', A = 5.0'). 2FGL J1022.7-5741, 
which is in the Westerlund 2 field, lies only 10' from 2FGL J1023.5— 5749c, so possibly 
the determination of its localization and/or localization uncertainty has been affected 
by this nearby source. 

PSR J1357-6429 near 2FGL J1356.0-6436 (9 99 , 9 = 9.1', A = 9.5'). 2FGL J1356.0-6436 
is a relatively isol ated source, but we note a p ossible association with the PWN 
HESS J1356-645 iLemoine-Goumard et al.lboilh . 



• PSR J1747-2958 near 2FGL J1747.1-3000 (6 99 . 9 = 2.9', A = 3.2'). 2FGL J1747.1-3000 
is located near the Galactic Center, and the localization of the source may be affected 
by systematic uncertainties in the diffuse Galactic emission model. 

In addition to the identified pulsars, four 2FGL sources are associated with radio pulsars: 

» 2FGL J1112.5-6105: PSR J1112— 6103 
» 2FGL J1632.4-4820c: PSR J1632-4818 
i 2FGL J171 7.5-5802: PSR J1717-5800 (?) 
» 2FGL J1928.8+1740c: PSR J1928+1746 (?) 



PSR J1717— 5800 has E = 2.3 x 10 32 , ten times lower than for any known 7-ray pulsar. The 
other three have E > 10 34 erg s -1 and the LAT team phase-folds 7 rays from their positions 
using radio rotation ephemerides as described by lSmith et al.l (120081 ) . Gamma-ray pulsations 
have not been detected for these pulsars. We mark two of the associations as questionable (?) 
because the corresponding 2FGL sources have spectra that are considerab ly softer (spectral 
index ~ 2.5) than typically observed for 7-ray pulsars ( lAbdo et al.ll2010tl ). 



The automatic association procedure also finds 21 2FGL sources to be associated with 
MSPs. Nineteen of those have unassociated counterparts in the 1FGL ca talog, and have been 



discovered in radio 1 


pulsar searches of unassociated 1FGL s 


Coenard et al. 


2011 


: 


Keith et al. 


2011; 


Hessels et al. 


2011) 



(e.g. 



Ransom et al. 2011 



enough to allow phase folding 7-rays from the directions of the newly discovered radio pulsars 
can require a year of radio observations to disentangle, e.g., binary orbital motion from annual 
parallax. As the ephemerides become available many of the unassociated 1FGL sources may 
reveal 7-ray pulsations, as has already occurred for several. Two 2FGL sources associated 
with MSPs have no 1FGL counterparts: 



• 2FGL J1023.6+0040: PSR J1023+0038 
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2FGL J1125.0— 5821: PSR J1125-5825 (Bates et al.ll201lh 



Tarn et al.l (120101 ) reported the LAT detection of 7-ray emission toward PSR J1023+0038, the 
only known rotation powered MSP in a quiescent LMXB. The spectrum of 2FGL J1023. 6+0040 
is rather soft (spectral index ~ 2.5) for an MSP, but the system i s sufficiently spec ial that 
this does not necessarily rule out the association (see discussion in lTam et al.ll2010l ). 



5.2.5. Pulsar wind nebulae 

Formally, we find 69 2FGL sources to be associated with PWNe, but except for three, 
all of them are also associated with known pulsars. Among those are three sources for 
which a dedicated analysis allowed us to identify both the pulsar and the PWN; they are 
summarized in Table [HJJ and the 2FGL catalog contains both the pulsar and the PWN as 
separate associated sources. For the other 63 2FGL sources, the observed pulsations firmly 
identify the pulsars as the primary source of the observed 7 rays, although some minor 
contribution from a PWN cannot be excluded. 

More interesting are the three PWN associations for which no pulsar has so far been 
identified. Those are: 

• 2FGL J1112.1-6040: G291. 0-0.1 

• 2FGL J1640.5-4633: G338.3-0.0 

• 2FGL J1745.6-2858: G359.98-0.05 (?) 

We mark the last association as questionable because this source is located in the imme- 
diate vicinity of the Galactic center where we know that the accuracy of our model of the 

Table 10. Identified PSR & PWN 





Pulsar 


Pulsar Wind Nebula 


Ref 




PSR 


2FGL 


PWN 


2FGL 




J0835- 


-4510 (Vela) 


J0835.3-4510 


Vela X 


J0833.1-4511e 


1 


J1509- 


-5850 


J1509.6-5850 


MSH 15-52 


J1514.0-5915e 


2 


J1826- 


-1256 


J1826.1-1256 


HESS J1825-137 


J1824.5- 1351c 


3 



References. - lAbdo et al ] feoioft : ^Abdo et al ] l|2010dh ; fcrondin et al. 



-80 - 



diffuse Galactic emission is intrinsically limited, and because the large density of potential 
counterparts makes a reliable source association difficult. 

In the 1FGL catalog we reported 6 sources associated with PWNs that were not also 
associated with known pulsars. Among those, two are among the 3 objects mentioned above 
(G338.3-0.0 and G359.98-0.05), one has turned out in fact to be a pulsar (2FGL J1135.3-6054), 
two are still unassociated 2FGL sources, but no longer associated with PWNs (2FGL J1552.8— 5609 
and 2FGL J1635.4— 4717c), and one no longer has a corresponding source in 2FGL (GO. 13— 0.11, 
see § Q . 



5.2.6. Globular clusters 



Eleven 2FGL sources are asso ciated with globu 



ar clusters. Among those, 9 have 
been published previously: 47 Tuc (lAbdo et al.ll2009ah . NGC 6266, N GC 6388, Terzan 5 , 



NGC 6440 NGC 6626, NG C 6652 flAbdo et all BOlOej ). Omega Cen (lAbdo et all EOlObh 
and M 80 ( ITam et al.ll2011bl ). In addition, we find two new associations: 



2FGL J1727. 1—0704: IC 1257. With an average significance of 4.1 this source is near 
the detection threshold. It is fitted using a power law with a spectral index of 2.2 ±0.1, 
yet a 3.5 o curvature significance may indicate that the spectrum is in fact curved. 

2FGL J1808.6-1950c: 2MS-GC01. This source has already been detected as 1FGL J1808.5-1954c, 
but the globular cluster catalog used for the association of 1FGL sources did not contain 
2MS-GC01, and consequently the source remained unassociated. 2FGL J1808.6— 1950c 
has an apparently-curved spectrum (3.9 cr significance) that is comparable to that of 
other globular clusters. 



Tarn et al.l (l2011bl ) have furthermore reported the detections of Liller 1, NGC 6139, 



NGC 6624, and NGC 6752 using LAT data. None of these clusters are formally associated 
with any of the 2FGL sources in the catalog. NGC 6624 is near 2FGL J1823.4-3014 (6 95 = 
7'7, A = 7'6), but the formal posterior association probability of 50% is below our adopted 
threshold. A source associated with NGC 6752 was in our initial list of seeds for the catalog; 
however, it did not pass the detection threshold of TS > 25 for the 2FGL catalog. We 
could not find evidence for any sources in our data that might be associated with Liller 1 or 
NGC 6139. 
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5.2.7. Supernova remnants 



SNRs are a special class in our association scheme because a substantial number of the 
known objects are sufficiently extended to be potentially resolved with the LAT. We thus use 
two separate strategies to search for SNR associations among the 2FGL sources. For SNRs 
with angular diameters < 20', i.e., SNRs that still should appear point-like to the LAT, 
we use the Bayesian scheme to search for associations. In total we find six 2FGL sources 
associated with point-like SNR, of which 2 are also associated with firmly identified pulsars. 
The remaining associations are: 



• 2FGL J1214.0-6237: G298.6-00.0 

• 2FGL J1911. 0+0905: G043.3-00.2 (aka W49B) 

• 2FGL J2022.8+3843c: G076.9+01.0 

• 2FGL J2323.4+5849: Gill. 7-02.1 (aka Cas A) 



None of them has a concurrent association with a PWN. Except for 2FGL J2022.8+3843c, 
all of them were already present and associated in 1FGL. 

In a second pass we search for all 2FGL sources for which the 95% confidence error 
radius overlaps with the (assumed) circular extension of the SNR. This provides a list of 
89 2FGL sources among which we estimate ~ 45% chance coincidences. Six of the 2FGL 
sources correspond to SNRs that were firmly identified as 7-ray sources based on their spatial 
extensions (IC 443, W28, W30, W44, W51C, and the Cygnus Loop), and 4 are the point- 
like SNRs listed above. Twenty of the 2FGL sources are firmly identified as being either a 
pulsar, a PWN, or a high-mass binary system. This leaves 59 2FGL sources that might be 
associated with an extended SNR, among which we expect ~ 26 chance coincidences. Due 
to this high chance coincidence rate, we do not claim any SNR association for this list of 
sources, but we give the 2FGL names and associations in Table [TT] for reference. 

Several of the SNRs have extensions that encompass multiple 2FGL sources (G132.7+01.3, 
Monoceros Loop, Pup A, Vela Junior, and G089.0+04.7), in which case the 2FGL sources 
might actually correspond to local maxima of extended emission regions. A number of the 
SNRs have been detected at TeV energies, which makes their possible detection also in the 
LAT energy range more plausible. Three 2FGL sources have concurrent PWN associations, 
which makes them also good pulsar or PWN candidates. We also note that one source, 
2FGL J2015. 6+3709, is likely to be variable, hence a physical association to CTB 87 is 
highly improbable. 
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In 1FGL, 4 1 7-ray sources we re listed in the corresponding table of overlaps with SNRs 
(see Table 7 of lAbdo et al.ll2010gl ). About half of SNRs that were found overlapping with 
1FGL sources are still in Table ITTj while the other half has not been found to overlap spatially 
with any of the 2FGL sources. This illustrates the relatively large uncertainty that is tied 
to these associations, and should present an additional warning to treat these potential 
associations with great care. 
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Table 11. Potential Associations for Sources Near SNRs 



2FGL name SNR name PWN name 



TcV name 



Common name 



J0128. 


.0+6330 


G127 


.1+00 


.5 


J0214. 


M 


-6251c 


G132 


.7+01 


.3 


J0218. 


.74 


-6208c 


G132 


.7+01 


.3 


J0221. 


.4-1 


-6257c 


G132 


.7+01 


.3 


J0503. 


.2+4643 


G160 


.9+02 


.6 


J0526. 


6+4308 


G166 


.0+04. 


.3 


J0538. 


.1+2718 


G180 


.0-01 


.7 


J0553. 


.9+3104 


G179 


.0+02 


.6 


J0631. 


.64 


-0640 


G205 


.5+00 


.5 


J0636. 


■ (H 


-0554 


G205 


.5+00 


.5 


J0637. 


8-t 


-0737 


G205. 


.5+00 


.5 


J0821. 


0- 


-4254 


G260. 


4-03 


4 


J0823. 


0- 


-4246 


G260 


.4-03 


4 


J0823. 


.4- 


-4305 


G260. 


.4-03 


.4 


J0842. 


9- 


-4721 


G263 


9-03 


3 


J0848. 


.5- 


-4535 


G266. 


2-01 


2 


J0851. 


.7- 


-4635 


G266. 


.2-01. 


.2 


J0853. 


.5- 


-4711 


G266. 


2-01 


.2 


J0855. 


4- 


-4625 


G266 


2-01 


.2 


J1112. 


.1- 


-6040 


G291. 


.0-00 


.1 


J1411. 


.9- 


-5744 


G315. 


.1+02. 


.7 


J1441. 


6- 


-5956 


G316 


3-00. 


.0 


J1521. 


.8- 


-5735 


G321 


.9-00. 


3 


J1552. 


.8- 


-5609 


G326. 


.3-01 


.8 


J1615 


0- 


-5051 


G332 


4+00. 


1 


J1628. 


.1- 


-4857c 


G335 


2+00. 


1 


J1631 


.7- 


-4720c 


G336 


.7+00. 


.5 


J1635. 


4- 


-4717c 


G337 


.2+00. 


.1 


J1640. 


.5- 


-4633 


G338. 


3-00 


.0 


J1712. 


.4- 


-3941 


G347. 


.3-00 


.5 


J1714. 


.5- 


-3829 


G348. 


5+00. 


1 


J1718. 


.1- 


-3725 


G350 


1-00. 


3 


J1727. 


3- 


-4611 


G343 


.0-06 


.0 


J1731. 


.6- 


-3234c 


G355 


.4+00. 


.7 


J1737. 


2- 


-3213 


G356 


3-00 


3 


J1738. 


.9- 


-2908 


G359 


.1+00. 


.9 


J1740. 


.4- 


-3054c 


G357 


.7-00. 


.1 


J1745. 


.5- 


-3028c 


G358. 


5-00. 


9 


J1745. 


.6- 


-2858 


G000. 


.0+00. 


.0 


J1802. 


3- 


-2445c 


G005. 


.4-01 


.2 


J1811. 


.1- 


-1905c 


G011. 


4-00. 


1 


J1828. 


3- 


1124c 


G020. 


.0-00 


2 


J1834. 


3- 


-0848 


G023. 


.3-00. 


.3 


J1834. 


.7- 


-0705c 


G024. 


7+00. 


6 


J1839. 


,7- 


-0334c 


G028. 


8+01. 


.5 


J1840. 


.3- 


-0413c 


G027. 


.8+00. 


.6 



RX J0852.0-4622 
RX J0852.0-4622 
RX J0852.0-4622 
RX J0852.0-4622 



G291. 0-0.1 



G338.3-0.0 



HESS J1616-508 



HESS J1634-472 
HESS J1640-465 
RX J1713.7-3946 
CTB 37A 



G359.98-0.05 



Monoccros Loop 
Monoccros Loop 
Monoceros Loop 
Pup A 
Pup A 
Pup A 
Vela 
Vela Junior 
Vela Junior 
Vela Junior 
Vela Junior 



Kes 25 
Kes 32 



HESS J1745-303 



HESS J1834-087 



CTB 37A 
RCW 114 

Tornado Nebula 

Sgr A East 
Bird 

W 41 
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8.4 GHz flux (mJy) 



Fig. 30. — Normalized histograms of the 8.4 GHz radio flux of CRATES sources (red: all 
sources, blue: objects associated with 2FGL sources). 
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5.2.8. Binaries 

The 2FGL catalog includes four high-mass binary systems, all of which have been firmly 
identified by their orbital modulation, and are described in separate publications: 



2FGL J0240.5+6113: LSI +61 303 flAbdo et al.ll2009ch . 



2FGL J1019.0-5856: 1FGL J1018.6-5856 flCorbet et al.ll201lh 



2FGL J1826.3-1450: LS 5039 flAbdo et al.ll2009d ). and 



2FGL J2032. 1+4049: Cygnus X-3 flAbdo et al.ll2009fl ). 



No further 2FG L source is associated with a high-mass X-ray binary from Liu's catalog 
( jLiu et al.l 120061 ) . All four sources were already present in the 1FGL catalog, yet the orbital 
modulatio n of 1FGL J1018.6— 5856 was only recently discovered in a blind search using the 
LAT data (lCorbet et al.lboilh . 



Formally, the automatic association procedure associates three 2FGL sources with low- 
mass X-ray binaries, but all three are located in globular clust ers, and the observe d emission 
can be readily explained by the combined emission of MSPs ( lAbdo et al.ll2009af ). We thus 
conclude that no low-mass X-ray binary systems have been identified in the LAT data after 
2 years of observations. We came to the same conclusion for 11 months of data in our study 
of the 1FGL associations ( lAbdo et al.ll2010gl ). 



5.2.9. Massive stars and open star clusters 

Among the massive star catalogs (O stars, Wolf-Rayet stars, Luminous Blue Variables) 
and the open cluster catalog we find only 2 possible associations with 2FGL sources: 



2FGL J1045. 0—5941: rj Carinae (LBV). The 7-ray emission of th i s well-known peculia r 
bin ary system has been studied in detail by lTavani et al.l (120091 ). lAbdo et al.l ( 12010hl ). 
and lFarnier et al.l (j201l[ ). yet a firm identification of the system through periodic orbital 
variability in 7 rays is still missing. 

2FGL J2030. 7+4417: HD 195592 (O star). This 09.51a type star is probably a short 
period (5.063 days) O+B binary s ystem at a distance of 1 .1 kpc that may have escaped 
from the open cluster NGC 6913 (IDe Becker et al.ll2010l ). We note, however, that the 
object is located in the Cygnus region where the high O star density easily could 
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lead to false associations and the complex diffuse emission may render precise source 
localization difficult. In addition, the spectral shape and the apparent lack of variability 
of 2FGL J2030. 7+4417 are similar to the characteristics of identified 7-ray pulsars. 
Hence we caution against overinterpreting this particular O star association and we do 
not list it in our final table. 



5.2.10. Multiwavelength associations 

In addition to the catalogs of classified sources, we also search for associations with 
catalogs of radio and TeV sources. Our association procedure for AGN heavily relies on 
associations with radio sources as most of the 7-ray emitting AGN are bright sources of 
radio emission (see § I5.2.2p . In fact, essentially all of the radio associations we find have 
been classified subsequently as AGN. 

Eighteen 2FGL sources that have not been associated with any object in one of our 
catalogs of known or plausible 7-ray-emitting source classes (our type 1 catalogs in § 15. 2j) 
have associations with extended TeV sources. However, due to the relatively large extents of 
the sources in the extended TeV catalog, we expect on average 20 false associations (cf. Table 
IH]), so from a statistical point of view, all 18 associations could be spurious. We discuss 2FGL 
associations with TeV sources more deeply in §5.31 



5.2.11. Other GeV Detections 



The automated association process compares the 2FGL source locations with other 
catalogs of sources seen at GeV energies. Results are shown in the main table for individual 
sources. From the Bright Gamma- Ray Source List (lAbdo et al.l l2009d[ ) we find 185 out of 
205 sources associated with 2FGL sources. Comparison with the 1FGL catalog was described 
in detail in § 14.21 In total we find 1099 out of 1451 1FGL sources that are associated with 
2FGL sources. 

The only contemporan eous catalog from a different instrument is the AGILE (1AGL) 
catalog ( jPittori et al.l 120091 ). which has 42 (out of 47) sources in common with the 2FGL 
catalog. The five 1AGL sources that are not formally associated (1AGL J0657+4554, 
1AGL J0714+3340, 1AGL J1022-5822, 1AGL J1803-2258 and 1AGL J1823-1454) all lie 
close to 2FGL sources and spatially overlap within their mutual 99% confidence localization 
uncertainties. Several 2FGL sources are associated with the same 1AGL source, and in total 
we find 57 2FGL sources associated with sources listed in 1AGL. 
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From the previous generation high-energ y 7-ray telescope, EG RET on the Compton 
Gamma Ray Observatory, the 3EG catalog (IHartman et al.l Il999f ) had 111 sources (out 
of 271) associated by the au tomatic process with 2FGL sources, while the EGR catalog 
( ICasandjian &: Grenierl 120081 ) had 66 (out of 188) sources associated with 2FGL sources. 
Also here we find several 2FGL sources that are associated with the same EGRET source. 
In total, 116 2FGL sources are associated with sources in 3EG, while 69 2FGL sources are 
associated with sources in EGR. The fractions of 3EG and EGR sources with 1FGL sources 



were s imilarly low and the discussion in the 1FGL catalog paper is still rel evant (lAbdo et al. 



2010g); we also refer the reader to a study of unassociated 1FGL sources (lAckermann et al. 



2012a). An EGRET catalog based on analysis of energies above 1 GeV (ILamb &: Macomb 



19971 ) found 46 high-confidence sources, of which 40 have clear 2FGL counterparts, 5 have 
close 2FGL sources just outside the 95% confidence contours, and only one (GEV 2026+4124 
in the confused Cygnus region) lacks a plausible 2FGL match. 

Through 2011 June, 94 flaring Fermi-LAT sources were detected and promptly re- 
ported in more than 150 Astronomer's Telegrams. Of these, 8 are not in 2FGL. For 6 
of these the flaring state was det ected outside the time interval covered by 2FGL: SBS 
0846 +513 (a new NLSyl sy stem: lDo~nato fc Perkir j |2 01l]) , SHBL J0013 55.9-185406 (see 



5.3: ISanchez fc Fegankoioh. PSR B1259-63 (see § 5J$; 



D'Ammando 



Abdo et al.ll2010a). PMN J1123-6417 



201 lh . PMN J1 913-3630 dDonato & Cheung! I2OI0L and the flar- 



(see § 4.2; 

ing source in the Galactic center region (jVasileiou et al.l 1201 ll ). The other two sources are: 
J1057-6027, flYasuda et aDl2009h detected in 2009 June, is not included in 1FGL and does 
not have a 2FGL counterpart but could be associa ted with 2FGL J1056.2 — 6021 using the 
99.9% confidence error radius; and PKS 1915—458 (jSokolovsky et al.ll2010l ) a faint and high 
redshift blazar [z = 2.47), detected in 2010 June, whose average flux between 2008 Au- 
gust and 2010 August is below the 2FGL catalog significance threshold. Also, we note that 
two 1FGL uniden tified flaring sources detected along the Galactic plane, 3EG J0903— 3531 
jHavs et alJliooj ) and J0910-5041 (2FGL J0910.4-5050 or 1FGL J0910.4-5055: lcheung etal 
2008|) are now associated with two unclassified AGN in 2FGL, PMN J0904-3514 an d AT20G 
J091 0— 5048 respectively. Furthermore, the 2FGL counterpart for J1512— 3221 ( [Wallace 
20101 ). which had no clear association, is 2FGL J1513.6— 3233 which is associated with blazar 
CRATES J1513-3234. 



5.3. TeV Source Associations 



2FGL sources that are positionally associated with sources seen by the ground-based 
TeV telescopes are of particular interest because the TeV band overlaps with the LAT energy 
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range, suggesting the potential for common emission mechanisms if the spectra match. As 
described in Tabled we investigated associations with the sources in the TeVCat compilation 
of detections. The compilation is growing with time, and information about the sources is 
subject to updates and refinements, but at any given time TeVCat represents a snapshot of 
current knowledge of the TeV sky. 

The association analysis was done separately for extended and point-like TeV sources, 
taking into account the statistical and systematic uncertainties in the source localization. 
The 'TeV column of Table H] lists associations with extended sources as 'E' and point-like 
sources as £ P\ As the table indicates, 85 2FGL sources are positionally consistent with 
TeVCat sources, although multiple 2FGL associations are seen for some TeV sources. In the 
FITS version of the catalog, we also provide the names of the associated TeV sources. 

Of the TeV sources considered for the associations performed here, most correspond 
to known objects at other wavelengths, in particular those that lie far from the Galactic 
plane. A large fraction (~ 50%) of the TeV Galactic sources, however, are still unidentified. 
Many of these have plausible counterparts while others remain unassociated despite deep 
searches for counterparts at other wavelengths. Among the firm identifications in the TeV 
regime, there are seven different source classes, and members of each of these source classes 
have been associated with 2FGL sources. In total, 85 TeV sources have 2FGL counter- 
parts. Eight of these TeV sources have more than one 2FGL association. RXJ0852. 0—4622 
has four 2FGL as sociations, and the following TeV sources have two each: Westerlund 1 



(lOhm . et all 120101 ) . Westerlund2, HESS J1632-478, RX J1713.7-3946, W28 f lAbdo et al. 



2010kh . HESS J1841-055, and MGRO J2019+37. One LAT source, 2FGL J2229. 0+6114 is 



associated with two TeV sources, Boomerang and G106.3+2.7. The LAT emission from two 
of the TeV sources, IC443 and MSH15— 52, is measured to be extended. 

The TeV class that has the most numerous associations with the 2FGL sources is the 



AGN class (see lAckermann et al.ll2011bl . for a more detailed discussion of the LAT AGN). 
There are currently 45 AGN detected at TeV energies and all but six of these are associated 
with 2FGL sources. The six that do not have 2FGL counterparts (SHBL J001355.9- 185406, 
1ES 0229+200, 1ES 0347-121, PKS 0548-322, 1ES 1312-423, and HESS J 1943+2130) are 
all high-frequency peaked BL Lacs. This is the subclass of AGN that tend to have the 
lowest bolometric luminosities and their second emission peaks at the highest energies. The 
six TeV AGN that did not reach the detection threshold to be included in the 2FGL catalog 
are among the weakest extragalactic TeV sources detected to date, ranging in flux from 



13 Thi s source has not been confir med to be a HBL but all available observations favor its classification as 
a HBL |HESS CollaboratioijEullbl) 
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0.4% — 2% of the flux of Crab Nebula at those energies. 



Both of t he starburst galaxi es detected at TeV energies, M82 and NGC253, have 2FGL 
counterparts (lAbdo et al.ll2010d ). 



Four high-mass binaries (HMBs) have confirmed detections in the TeV regime. Two 
of these, LSI +61 303 and LS 5039, hav e 2FGL counterparts and are already the subject 
of LAT publications fUbdo et al.ll2009cl l3). We note that, although not in 2FGL, the TeV 
binary PSRB1259-63 has been detected by the LAT jTam et al.lboilal : lAbdo et alJboilah . 
This system is a radio pulsar in orbit around a Be star with an orbital period of ~ 3.4 
years. During the time span of the 2FGL data, the system was far from periastron and no 
significant GeV emission was detected, but when the system approached periastron, variable 
emission, including flaring behavior, was observed by the LAT. 

The PWNe comprise the second most numerous identified TeV class that is associated 
with 2FGL sources; of the 25 PWNe in TeVCat, 16 are associated with 2FGL sources. 
Indeed, the association betwe en GeV 7-ray PSRs and t he PWNe visible in the regime has 
been well established already ( lAbdo et al.ll2010tl 12012a ). 



During the second year of LAT data taking, many more superno ya remnants (SNRs ) 
known at TeV ene rgies were detected at GeV ener gies such as Cas A (lAbdo et al.ll2010ml ). 
RX J1713. 7-3946 fUbdo et alihoildh . and Vela Jr franaka et alkoilh . Of the five SNR/Molecular 
Cloud associations in TeVCat, all but one (G318. 2+0.1) have been associated with 2FGL 
sources. Ten shell-type SNR have been detected at TeV energies and five of these now have 
2FGL counterparts so, the GeV- TeV association is established although there are still many 
open questions. Fermis non-detection of RCW 86 is surpri sing since it is one of th e brightest 
TeV SNR, with a flux of ~ 10% that of the Crab Nebula jAharonian et al.ll200gh . 



Sources of particular interest are those that are positionally consistent between the LAT 
and TeV telescopes but have no obvious associations with objects at longer wavelengths. 
Among the TeV sources that have no clear identifications, 17 are associated with 2FGL 
sources. In addition to these, although not formally associated with LAT sources using 
the automatic pipeline (§ I5.2|) . some other TeV sources have possible 2FGL counterparts, 
for example, HESS J1843— 033, which has two potential 2FGL counterparts. Establishing a 
physical connection through spectral or variability studies may help determine the nature of 
these sources. In many cases, a GeV counterpart could prove crucial for our understanding 
of the nature of the TeV source, in particular for the following objects: 

• 2FGL J1022.7-5741 and 2FGL J1023. 5-5749 are spatially consistent with HESS J1023- 
575, itself not yet firmly identified, but noted for its possible connection to t he young 



stellar cluster Westerlund 2 in the star-forming region RCW 49, as discussed by lHESS Collaboration 
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bond). 



2FGL J1427.6-604 8 is associated with HES S J1427-608 which is, so far, without plau- 
sible counterparts (jAharonian et al.ll2008bl ). 



2FGL J1503.9— 5800 is spatially coincident with the TeV source HESS J1503-582, 
which is tentatively associated with a forbidden velocity region of interstellar gas 
jRenaud et al.ll2008h . 



2FGL J1507. 0—6223 is spatially consistent with HESS J1507-622 flHESS Collaboration 
2011al ). so far the only TeV unidentified source that is markedly offset from the Galactic 
plane (~ 3?5). 

2FGL J1615.2— 5138 is spatially consistent with one of the brightest (~ 25% of the Crab 



Nebula flux) TeV unidentified sources, HESS J1614— 518 (jAharonian et al.ll2006l ) 



2FGL J1650.6— 4603 is spatially associated with a TeV source tentatively associated 
with the Westerlund 1 star-forming region. (The other 2FGL source that is spa- 
tially associated with this TeV source is 2FGL J1648.4— 4612, which is pulsar PSR 
J1648-4611.) 

2FGL J1848.2— 0139: this source is consistent with the TeV source, HES S J1848-018, 



which is suspected to be correlated with the star-forming region W43 (jChaves et al. 
2008T ). 



As discussed in § 12.21 the Galactic Center region is particularly complex and its study is 
beyond the purpose of this paper; we do, however, find possible associations with all of the 
TeV 7-ray sources detected in this region, although not all were formally a ssociated by the au- 
tomatic pipeline analy sis: the Galactic Center source (Ucero et alJboogh. HESS J1745 -303 



(2FGL J1745.5-3028c; lAharonian et al.ll2008ah. HESS J1741-302 



HESS J1747-248 (2FGL J1748.0-2447; IHESS Collaborationll2011df ) 



Tibolla et al.ll2008f ) and 



5.4. Properties of Unassociated Sources 

Among the 1873 sources in the 2FGL catalog, 575 (31%) remain unassociated. Their 
distribution on the sky is compared in Figure I3T1 to the distribution of the associated sources. 
We note a number of interesting features in the map that should be kept in mind when 
considering unassociated 2FGL sources. 
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Table 11 


— CJontinupd 


2FGL name 


SNR name PWN 


name TeV name Common name 


J1841.2-0459c 


G027.4+00.0 


Kes 73 


J1849.3-0055 


G031. 9+00.0 


Kes 77, 3C 391 


J looU. I — UU14C 






J1852.7+0047c 


G033. 6+00.1 


Kes 79 


J1916.1+1106 


G045.7-00.4 




J1932.1+1913 


G054.4-00.3 




J2015.6+3709 a 


G074.9+01.2 


CTB 87 


J2019. 1+4040 


G078. 2+02.1 


VER J2019+407 Gamma Cygni 


J2041. 5+5003 


G089. 0+04.7 




J2043.3+5105 


G089.0+04.7 




J2046.0+4954 


G089. 0+04.7 




J2333.3+6237 


G114.3+00.3 




J2358.9+6325 


G116. 5+01.1 





a Source is likely to be variable. 




Fig. 31. — Sky distribution of associated (dots) and unassociated sources (large symbols). 
Sources that were nagged are marked by a plus. In particular, we mark variable unassociated 
sources (TS var > 41.6) using filled blue circles, unassociated sources with a curved spectrum 
(Signif _Curve > 4) by red crosses, and all other unassociated sources by open black circles. 
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First, the number of unassociated sources decreases with increasing Galactic latitude. 
This is best illustrated by a latitude histogram of the fraction of unassociated 2FGL sources, 
shown in Figure [32J We plot here the data as function of the sine of Galactic latitude as in 
this representation an isotropic distribution will appear as a flat profile. In contrast to that, 
we find that the fraction of unassociated sources decreases with latitude, with the decrease 
being steeper at positive latitudes. This asymmetry is also present in the absolute numbers: 
above Galactic latitudes b > 60° only 3 sources in 2FGL are unassociated, while below 
b < —60° we find 12 unassociated sources. This may be due to relative completeness in 
the north vs. the source of the counterpart catalogs used for the source association analysis 

(§E2D. 

Second, the numbers of unassociated sources increase sharply below |6| w 10°. This 
is attributable to the relative lack of sources below |6| < 10° in many of the extragalactic 
source catalogs that we use for source association. The Milky Way is a bright source of radio 
emission, limiting sensitive searches for extragalactic sources near the Galactic plane. Fur- 
thermore, optical identifications of radio sources are hampered by the important interstellar 
obscuration, leaving many radio sources unclassified. 

Third, the numbers of 2FGL sources with curved spectra increase at low Galactic lat- 
itudes, as can be seen in the latitude histogram (dashed line in Fig. [32]) and the sky map 
(red crosses in Fig. |3T]) . The sky map indicates that these sources tend to cluster in regions 
of bright Galactic diffuse emission, such as the inner Galactic ridge (Galactic longitudes 
330° < / < 30°), the Cygnus region (I « 80°), the Norma spiral arm tangent (I « 330°) or 
the Crux spiral arm tangent (/ ~ 300°). Whether this clustering is diagnostic of the physical 
natures of the sources, or whether it indicates systematic uncertainties in the Galactic diffuse 
emission model that resulted in spurious source detections remains a possibility. We note, 
however, that the fraction of sources with curved spectra among the unassociated sources 
is greater (28%) than the fraction of curved spectra sources among the associated sources 
(16%). Because the spectrum of the Galactic diffuse emission at low latitudes is itself well 
represented with a curved spectrum, at least some fraction of the unassociated 2FGL sources 
at low latitudes may be local emission maxima of diffuse Galactic emission that are not ade- 
quately modeled by our Galactic diffuse model; see the discussion in § 13.91 and the definitions 
of the several analysis flags that are related to the model of the Galactic diffuse emission in 

Fourth, a substantial fraction of the unassociated sources have at least one analysis flag 
( §3.10p set. We find that 51% of the unassociated sources have been flagged due to various 
issues, while only 14% of the associated sources have been flagged. None of the flags is 
related to our association procedure itself, but they identify a number of conditions that can 
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shed doubt on the physical reality or localization quality of a source. The fact that such a 
large fraction of unassociated sources are flagged may indicate that some of these sources 
are indeed not real. We emphasize that the analysis flags should be taken into consideration 
when using the 2FGL catalog. 

Fifth, 25 unassociated sources (4%) have been flagged as variable, and the spatial distri- 
bution of these sources appears rather isotropic. These sources are good candidates for being 
as-yet unassociated AGN, as this is the source class that shows the largest flux variability 
in LAT data. 



6. Conclusions 

The second Fermi LAT catalog is the product of a comprehensive analysis of the first 
2 years of LAT science data. In several ways it is an advance over the 1FGL catalog, which 
was based on the first 11 months of data. The 2FGL analysis takes advantage of the new 
P7_V6 Source event selection and IRFs, which in particular provide increased effective area 
in the range below ~200 MeV. The analysis also uses a refined model for the Galactic diffuse 
emission. The source detection and localization analyses were advanced for the 2FGL anal- 
ysis to iteratively optimize the definitions of the 'seed' sources used for the final likelihood 
analysis step. Both analysis steps allowed for non-power-law source spectra and also incorpo- 
rated special models for spatially extended sources. The source association analysis was also 
extensively updated for the 2FGL catalog, with updated catalogs of counterparts and local 
determinations of counterpart densities. For AGN, the association analyses also included 
methods that took into account radio and X-ray properties of potential counterparts. 

The 2FGL catalog contains 1873 sources. In developing the catalog analysis, we re- 
evaluated a number of the analysis flags used to tag sources with unusual or potentially 
problematic properties. The most-prominent flag is the 'c' designator, which we have ap- 
pended to the names of 162 sources, that indicates potential confusion with interstellar 
diffuse emission or an artifact in the model for the diffuse emission. A number of other flags 
are defined, and 315 sources have one or more of these other flags set. 

The 2FGL catalog represents a new milestone in high-energy 7-ray astrophysics. As 
with any astronomical catalog, 2FGL enables a wide range of astrophysical research. For 
individual objects, the spectra and light curves offer opportunities for multiwavelength mod- 
eling that can lead to better physical understanding of sources. The catalog as a collection 
allows population studies for 7-ray-only sources and for comparative studies with other 
wavelengths. In the catalog, 127 sources are considered to be identified, and plausible asso- 
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ciations are proposed for more than 1000 AGN. In all identifications or associations of 2FGL 
sources with 15 classes of counterparts are proposed. In addition the fact that 575 of the 
2FGL sources have no plausible counterparts among known 7-ray-producing source classes 
presents discovery opportunities similar to those already found with the Fermi LAT Bright 
Source List and 1FGL catalog. Even the absence of 2FGL sources in predicted source classes 
such as clusters of galaxies will stimulate additional research into why these known sources 
of nonthermal radiation are not producing 7 rays at a level yet detectable with the LAT. We 
look forward to extensive use of this catalog in high-energy astrophysics. 



We dedicate this paper to the memory of our colleague Patrick Nolan, who died on 
2011 November 6. His career spanned much of the history of high-energy astronomy from 
space and his work on the Large Area Telescope began nearly 20 years ago when it was just 
a concept. Pat was a central member in the operation of the LAT collaboration and he is 
greatly missed. 

The Fermi-LAT Collaboration acknowledges generous ongoing support from a number 
of agencies and institutes that have supported both the development and the operation of the 
LAT as well as scientific data analysis. These include the National Aeronautics and Space 
Administration and the Department of Energy in the United States, the Commissariat a 
l'Energie Atomique and the Centre National de la Recherche Scientifique / Institut National 
de Physique Nucleaire et de Physique des Particules in France, the Agenzia Spaziale Italiana 
and the Istituto Nazionale di Fisica Nucleare in Italy, the Ministry of Education, Culture, 
Sports, Science and Technology (MEXT), High Energy Accelerator Research Organization 
(KEK) and Japan Aerospace Exploration Agency (JAXA) in Japan, and the K. A. Wallen- 
berg Foundation, the Swedish Research Council and the Swedish National Space Board in 
Sweden. 

Additional support for science analysis during the operations phase is gratefully acknowl- 
edged from the Istituto Nazionale di Astrofisica in Italy and the Centre National d'Etudes 
Spatiales in France. 

This work made extensive use of the ATNF pulsar catalog^ ( Manchester et al. 2005 ). 



This research has made use of the NASA/IPAC Extragalactic Database (NED) which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. 
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A. Localization Power 



In the maximum likelihood formalism, the error ellipse (§ 13. 1.4j) is given by the covari- 
ance matrix of the position parameters after the fit. One can obtain an approximate but 
reasonably accurate estimate of the localization power of the Fermi-LAT for a point source, 
assuming that the diffuse background is locally uniform and considering only one source. In 
that approximation the error ellipse is a circle and the 1 a localization precision of a source 
along any direction A#o is given by 



d 2 log C 


i 


d 2 TS 


del 


~ 2 


del 



(Al) 



and is relate d to the 95% error radius by r 95 /A^ = a/— 2 log(0.05) = 2.45. Along the lines 
of Eq. Al of lAbdo et al.l (j2010gh . denoting S(E) the source spectrum, B(E) the background 
spectrum per unit solid angle, T the equivalent on-axis observing time, A e fi(E) the on-axis 
effective area and the local source to background ratio g(e,E) = S(E)PSF(e, E)/B(E) one 
may write 



A# - 2 = T 



Wi(E)dlogE 



Wi(E) = nEA eS {E)- 



S(E) 2 r /<9PSF\ 2 sin 9d6 



B(E) Jo \ de J l + g(9,E) 



(A2) 
(A3) 



after integrating the cos 2 <fi term arising from the projection along one direction. Here Wi(E) 
is the contribution to A6*q 2 per unit log(energy). It is illustrated in Figure 1331 for a power- 
law source spectrum at high latitude. Not surprisingly, the localization depends even more 
on high energy (w here the core PSF is narrowest) than the detection itself (Figure 18 of 
Abdo et al.ll2010gl ). For that reason, the average effect of confusion on localization is small, 
because it is important only at those energies when the average angular distance between 
sources (2?8 at high latitude) is comparable to the PSF width. 

For each spectral index it is possible to compute the detection threshold and then the 
localization precision at the detection threshold. This is normally the worst error radius 
one may expect in the catalog. That prediction is compared on Figure [34] with the actual 



This preprint was prepared with the A AS IATgX macros v5.2. 
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95% error radius. The curve accurately predicts the dependence on spectral index. The 
localization is worse for softer sources, but only by a relatively small factor at a given TS. 
A few sources are above the line. This can happen for purely statistical reasons, because 
the background and exposure depend a little on direction even after taking out the Galactic 
plane, or because of another nearby source. The highest point (worst error ellipse) is 2FGL 
J1952.6— 3252 which is specifically flagged for imperfect localization (Flags 8 and 9 set). 



B. Quality of the Fit 

In order to illustrate the global quality of the main spectral fit (§ 13. 2p . we show in 
Figures [35] and [36] the spatial and spectral residuals over a large sky region rather than an 
individual Rol which could hide cross-talk issues. We chose the Galactic anticenter which is 
halfway between the quiet high-latitude regions and the most difficult Galactic Ridge regions 
discussed in § 13.91 

We fit the same parameters as in an ordinary Rol: normalizations of the isotropic and 
Galactic components Ki SO and K ga i, and corrective slope of the Galactic component T ga i, 
such that the correction to the Galactic model is K ga i(E/E )~ r 9 al with E set to 500 MeV. 
The fitted parameters were i\~ iso = 0.973, K ga i = 1.003 and T ga i = 0.029. 

The spatial residuals are scaled to the Poisson noise in each pixel in order to quantify 
whether the deviations are significant. What is shown is (data — model) / V model. The 
pixel size is large enough that there are about 100 counts per pixel outside the plane. The 
distribution of spatial residuals on 0?5 pixels follows very closely a normal law. Its standard 
deviation is only 1.1 a, implying that the intrinsic fluctuations are about 0.5 a, or 5%. They 
appear to be on a scale of a few degrees. The spectral residuals are a few percent and 
evolve slowly with energy. Those are small imperfections of the diffuse model, which show 
up because of the very high statistical quality of the data. Their impact on sources is limited 
because the residuals are on a larger scale than the LAT PSF except at low energy. It is 
quantified in § 13.71 



C. Description of the FITS Version of the 2FGL Catalog 

The FITS format version of the 2FGL catalog^] has four binary table extensions. The ex- 
tension LAT_Point_Source_Catalog Extension has all of the information about the sources, 



3 The file is available from the Fermi Science Support Center, http://fermi.gsfc.nasa.gov/ssc 
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including the monthly light curves (Tab. [T2|) . 

The extension Hist_Start lists the Mission Elapsed Time (seconds since 00:00 UTC on 
2000 January 1) of the start of each bin of the monthly light curves. The final entry is the 
ending time of the last bin. 

The extension GTI is a standard Good-Time Interval listing the precise time intervals 
(start and stop in MET) included in the data analysis. The number of intervals is fairly 
large because on most orbits (~95 min) Fermi passes through the South Atlantic Anomaly 
(SAA), and science data taking is stopped during these times. In addition, data taking is 
briefly interrupted on each non-SAA-crossing orbit, as Fermi crosses the ascending node. 
Filtering of time intervals with large rocking angles, other data gaps, or operation in non- 
standard configurations introduces some more entries. The GTI is provided for reference 
and would be useful, e.g., for reconstructing the precise data set that was used for the 1FGL 
analysis. 

The extension ExtendedSources contains information about the 12 spatially extended 
sources that are modeled in the 2FGL catalog, including locations and shapes (Tab. [T31) . 
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Table 12. LAT 2FGL FITS format: LAT_Point_Source_Catalog Extension 



Column 


Format 




Unit 


Description 


Source_Namc 


18A 








RAJ2000 


E 


deg 




Right Ascension 


DEJ2000 


E 


dcg 




Declination 


GLON 


E 


dcg 




Galactic Longitude 


GLAT 


E 


dcg 




Galactic Latitude 


Conf_68_ScmiMajor 


E 


dcg 




Long radius of error ellipse at 68% confidence 


Conf_68.SemiMinor 


E 


dcg 




Short radius of error ellipse at 68% confidence 


Conf_68_PosAng 


E 


dcg 




Position angle of the 68% long axis from celestial North, 
positive toward increasing RA (eastward) 


Conf_95_SemiMajor 


E 


dcg 




Long radius of error ellipse at 95% confidence 


Conf_95_SemiMinor 


E 


dcg 




Short radius of error ellipse at 95% confidence 


Conf_95.PosAng 


E 


dcg 




Position angle of the 95% long axis from celestial North, 
positive toward increasing RA (eastward) 


Signif_Avg 


E 






Source significance in a units (derived from Test Statistic) 


Pivot_Encrgy 


E 


MeV 




Energy at which error on differential flux is minimal 


Flux_Dcnsity 


E 


cm 


MeV" 1 s- 1 


Differential flux at Pivot_Energy 


Unc_Flux_Dcnsity 


E 


cm -2 


MeV" 1 s- 1 


1 a error on differential flux at Pivot_Encrgy 


SpectraLIndex 


E 






Best fit photon number power-law index. For LogParabola spectra, 
index at Pivot_Energy; for PLExpCutoff spectra, low energy index. 


Unc_Spcctral_Indcx 


E 






1 a error on SpectraLIndex 


FluxlOOO 


E 


cm -2 


s- 1 


Integral flux from 1 to 100 GeV 


UncFluxlOOO 


E 


cm 


s- 1 


1 a error on integral flux from 1 to 100 GeV 


Energy .FluxlOO 


E 


erg cm 2 s 1 


Energy flux from 100 MeV to 100 GeV obtained by spectral fitting 


Unc.Energy_FluxlOO 


E 


erg cm -2 s _1 


1 cr error on energy flux from 100 MeV to 100 GeV 


Signif.Curve 


E 






Significance (in a units) of the fit improvement between power-law 

and either LogParabola (for ordinary sources) or PLExpCutoff (for pulsars). 

A value greater than 4 indicates significant curvature. 


SpcctrumType 


18A 






Spectral type (PowerLaw, LogParabola, PLExpCutoff). 


beta 


E 






Curvature parameter (/3) for LogParabola. NULL for other spectral types 


Unc_beta 


E 






1 <x error on /3 for LogParabola. NULL for other spectral types 


Cutoff 


E 


MeV 




Cutoff energy as exp(-E/Cutoff) for PLExpCutoff. NULL for other spectral types 


Unc.Cutoff 


E 


MeV 




1 a error on cutoff energy for PLExpCutoff. NULL for other spectral types 


PowcrLawJndex 


E 






Best fit power-law index. Equal to SpectraLIndex if SpecrumType is PowerLaw. 


Flux30_100 


E 


cm~ 2 


s- 1 


Integral flux from 30 to 100 MeV (not filled) 


Unc_Flux30.100 


E 


cm~ 2 


s- 1 


1 cr error on integral flux from 30 to 100 MeV (not filled) 


Sqrt.TS30.100 


E 






Square root of the Test Statistic between 30 and 100 MeV (not filled) 


Fluxl00.300 


E 


cm~ 2 


s- 1 


Integral flux from 100 to 300 MeV 


Unc.Fluxl00.300 


E 


cm -2 


s- 1 


1 a error on integral flux from 100 to 300 MeV a 


Sqrt.TS100.300 


E 






Square root of the Test Statistic between 100 and 300 MeV 


Flux300.1000 


E 


cm~ 2 


s- 1 


Integral flux from 300 MeV to 1 GeV 


Unc.Flux300.1000 


E 


cm 


s- 1 


1 cr error on integral flux from 300 MeV to 1 GeV a 


Sqrt.TS300.1000 


E 






Square root of the Test Statistic between 300 MeV and 1 GeV 


Fluxl000.3000 


E 


cm~ 2 


s- 1 


Integral flux from 1 to 3 GeV 


Unc.Fluxl000.3000 


E 


cm~ 2 


s- 1 


1 a error on integral flux from 1 to 3 GeV a 


Sqrt.TS1000.3000 


E 






Square root of the Test Statistic between 1 and 3 GeV 


Flux3000.10000 


E 


cm~ 2 


s- 1 


Integral flux from 3 to 10 GeV 
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32. — Latitude distribution of unassociated sources. 
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Fig. 33. — Theoretical contribution (Wi(E) of Eq. IA3I) to A9 2 per Ms and per \og(E) 
interval as a function of energy for a TS = 100 power-law source over the average background 
at \b\ > 10°. The assumed photon spectral index is 2.2. 



- 107 - 




Fig. 34. — 95% error radius of sources at \b\ > 10° as a function of spectral index. The 
line shows the theoretical error radius for an isolated source at the detection threshold of 
TS = 25 over the average extragalactic background. 
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-5-3-11 3 5 



Fig. 35. — Residuals (in a units) in 0?5 pixels over a 60 x 60° area around the Galactic 
anticenter, summed over the full energy range (100 MeV to 100 GeV). All sources were fixed 
to the catalog values and the diffuse parameters were fitted as in an ordinary Rol (§ 13. 2(1 . 
The pixels used in the source fitting process were much smaller. The larger pixels used here 
allow reducing the statistical fluctuations to 5% in the Galactic plane and 10% at the top 
and bottom of the plot. 
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Fig. 36. — Left: Fit to the full spectrum integrated over the same anticenter region as in 
Figure 1351 The spectral bins are the same as in the source fitting process. The dotted, dashed 
and dash-dotted lines are the Earth limb, isotropic and Galactic components, respectively. 
The asterisks show the total source contribution (dominated by the Geminga and Crab 
pulsars). The full line is the sum of all model contributions, to be compared with the data 
(plus signs). The statistical errors on the data are shown but barely visible except at high 
energy. Right: Fractional residuals (data/model— 1) with statistical error bars. The residuals 
are statistically significant because of the very large number of events (2.8 x 10 6 over that 
area) but are only a few percent. 
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Table 12— Continued 



Column 


Format 


Unit 


Description 


Unc_Flux300CL10000 


E 


cm 


2 s" 1 


1 a error on integral flux from 3 to 10 GeV a 


oc|rL_± oouuu-iuuuu 


Hi 






Scjuarc root of the Test Statistic between 3 and 10 GcV 


r 1UX1UUUU LUUUUU 


IT 
£j 


cm" 


2 S" 1 


integral mix irorn iu 10 iuu \j& v 




J-T 


cm" 


2 s" 1 


1 (7 error on integral flux from 10 to 100 GeV a 


c nr + Tmnnnn innnnn 


E 






Scjuare root of the Test Statistic between 10 and 100 GcV 


V ouL IcxlJlll by _lilQCX 


E 






Sum of 2 x Log(Likelihood J comparison between the flux fitted in 24 time 
segments and a flat light curve over the full 2-ycar catalog interval. 

A \. r '"} 1 1 1 f> 1 VT ' +~ O 1 ' f~ r 1 *"1 t' 1 1 1 (~'\ 1 1 » 1 ( 1 1 i'"' '"1 +~ OO 1 ft- C M !^ 1 1 f'O /"^T llOl Tl ()' Q t" iH "\r Gl^l i vr*/ - * 
jTY. VrtlLlL, til LlleLIl t ±_l_.lJ-± 1I1L1.11^>£L \ _L /0 LTIOjIIUL. L>1 [JL-1I1H, d a LLdAJ. V 


Signif_Peak 


E 






Source significance in peak interval in a units 


Flux_Pcak 




cm" 


2 s" 1 


rcaK integral nux irom iuu ivicv ro iuu taev 


T Jy-i r - TVl 11V T^/~i \r 

U I1L._JT 1 LIX _i c t±K 




cm" 


2 s" 1 


1 (7 error on peak integral flux 


Time_Peak 


n 


s (MET) 


Time of center of interval in which peak flux was measured 


Peak_Intcrval 


FT 


s 




Length of interval in which peak flux was measured 


Flux_History 


111? 


cm" 


2 s" 1 


Integral flux from 100 MeV to 100 GeV in each interval (best fit from 
likelihood analysis with spectral shape fixed to that obtained over 2 years). 


Unc_Flux_History 


HE 


cm" 


2 s" 1 


Error on integral flux in each interval using method 

indicated in Unc_Flag_History column and added in quadrature 

with 3% systematic component. 


Unc_Flag_History 


11B 






1 if it is half of the difference between the 2 a upper limit 

and the maximum-likelihood value given in Flux_History, if it is the 

1 <t uncertainty derived from a significant detection in the interval 


Extended_Source_Name 


18A 






Cross-reference to the ExtendedSources extension for extended sources, if any 


0FGL_Name 


18A 






Name of corresponding 0FGL source, if any 


lFGL_Name 


18A 






Name of corresponding 1FGL source, if any 


ASSOC.GAM1 


18A 






Name of likely corresponding 1AGL source 


ASSOC_GAM2 


18A 






Name of likely corresponding 3EG source 


ASSOC_GAM3 


18A 






Name of likely corresponding EGR source 


TEVCAT_FLAG 


A 






P if positional association with non-extended source in TeVCat 

E if associated with a more extended source in TeVCat, N if no TeV association 


ASSOC_TEV 


24A 






Name of likely corresponding TeV source from TeVCat 


CLASS1 


3A 






Class designation for associated source; see Table [6] 


CLASS2 


3A 






Second class designation for associated source 


ASSOC1 


24A 






Name of identified or likely associated source 


ASSOC2 


24A 






Alternate name of identified or likely associated source 


Flags 


I 






Source flags (binary coding as in Table [3]l 



a The upper limit is set equal to if the flux in the corresponding energy band is an upper limit (TS < 10 in that band). The 
upper limits are 2 a. 
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Table 13. LAT 2FGL FITS format: ExtendedSources Extension 



Column 


Format 


Unit 


Description 


Source_Namc 


18A 






lFGL_Name 


18A 






RAJ2000 


E 


dcg 


Right Ascension of ccntroid 


DECJ2000 


E 


dcg 


Declination of ccntroid 


GLON 


E 


deg 


Galactic Longitude of centroid 


GLAT 


E 


dcg 


Galactic Latitude of centroid 


ModeLForm 


24A 




Spatial shape (2D Gaussian, Disk, Ring, Template, ...) 


ModeLSemiMajor 


E 


dcg 


Long radius of source. Full size for bounded shapes (disk, ring). 
68% containment for unbounded shapes (Gaussian) 


ModeLSemiMinor 


E 


dcg 


Short radius of source 


ModeLPosAng 


E 


dcg 


Position angle of the long axis from celestial North, 
positive toward increasing RA (eastward) 


SpatiaLFilcnamc 


68A 




Name of spatial template file a 



SpatiaLFilcname refers to external files that should be included with the catalog distribution. 



